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pi  ease  note  the  following  corrections: 
Page  25. 


Page  28. 
Page  44, 

Page  50. 
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Page  52 . 


Page  53. 


Fig.  15.  The  arrows  of  the  center  figure  should  be 
reversed,  and  the  phrase  (Advance  of  Line  of  Nodes) 
should  be  deleted. 

Footnote.  Change  Move  strictly  to  *More  strictly. 

First  equation.  (d,  |3 ,  0„,  y)  should  be  changed  to 
(a,  |3,  0g,  y)  . 

3(a).  Add  2i r  as  in  1(a). 
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ABSTRACT 


This  report  assembles  the  results  of  studies  in  the  field  of  space 
flight  thermal  radiation  simulation  to  provide  a  system  to  compute 
real-time  information  about  the  intensity  and  distribution  of  thermal 
radiation  fields  illuminating  a  vehicle  in  space.  The  details  are 
limited  to  vehicles  in  near  earth  orbits,  but  the  discussion  is  suffi¬ 
ciently  generalized  to  indicate  how  transfer  trajectories  and  alternate 
primary  bodies  may  be  incorporated.  The  time  functions  obtained 
may  eventually  control  the  sequencing  of  parameters  in  the  Mark  I 
aerospace  environment  simulator  or  at  least  provide  an  independent 
means  to  evaluate  system  performance  when  the  simulator  becomes 
operational. 


PUBLICATION  REVIEW 


This  report  has  been  reviewed  and  publication  is  approved. 


G.  M.  Arnold 

Chief,  Engineering  Division 
Space  Systems  Office 


,^/John  D. 

Lt  Col,  USAF 

Chief,  Space  Systems  Office 


111 


AEDC-TDR-63-206 


CONTENTS 

Page 

ABSTRACT . .  iii 

NOMENCLATURE .  vii 

1.  0  INTRODUCTION .  1 

2.0  THE  GENERAL  COMPUTER  PROBLEM . '  4 

3.0  SOLAR  CONTROL  AND  ORBITS 

3.  1  Orbit  Definition .  7 

3.  2  Orbit  Calculations .  8 

3.  3  Coordinate  System .  11 

3.  4  Orientation .  12 

3.  5  Roll  and  Pitch .  14 

3.  6  Orbits .  18 

3.  7  Solar  Bank  Control . .  .  27 

4.  0  SECONDARY  ILLUMINATION .  28 

4.  1  Secondary  Illumination  Simulation .  30 

4.  2  "Birdcage"  Patterns  and  Matching .  34 

4.  3  Illumination  Integrals .  40 

4.  4  Radiance  Functions .  46 

5.0  CONCLUSIONS . 47 

REFERENCES .  47 

APPENDIXES 

I.  A  List  of  7074  Programs  Concerned  with 

Mark  I  Illumination  Simulation .  49 

II.  Special  Cases  of  Albedo  Integral .  50 

III.  Earth  Radiance  Integral .  54 

IV.  Sample  Profile .  59 

TABLES 

1.  Program  Block  Outline .  5 

2.  Program  Elements .  7 

ILLUSTRATIONS 

Figure 

1.  The  Simulation  Problem .  1 

2.  Illumination  of  Satellite .  2 

3.  Typical  Shapes  and  Properties .  3 


v 


AEDC-TDR-63-206 


Figure  Page 

4.  Distinction  between  Rotational  Quantities  and 

Revolutional  Quantities  in  Flight .  8 

5.  Calculation  of  Node  Vectors .  9 

6.  Earth- Centered  Inertial  Coordinate  System  .  ...  11 

7.  Orientation  Types .  13 

8.  Roll  and  Pitch  Geometry .  15 

9.  Roll  and  Pitch  Functions  of  Orbit  Position  ....  17 

10.  Roll  Rate  Editing .  17 

11.  Eclipse  Geometry . 19 

12.  Eclipse  Test . 20 

13.  Orbit  Coordinate  Angles .  22 

14.  Orbit  in  Space  and  in  Plane  Illustrating  the 

Euler  Angles,  0.,  ,  « .  22 

15.  Secular  Perturbations  Caused  by  Earth 

Oblateness .  25 

16.  Geometry  for  Study  of  Longitudinal  Flux 

Distributions .  30 

17.  Longitudinal  Flux  Distributions  for  Line  Source 

Arrays  as  Function  of  Distance  from  Line  ....  32 

18.  Approximate  Effect  of  Reflecting  End  Baffles  on 
Longitudinal  Flux  Distribution  for  a  Single  Line 

Source  Element .  3  3 

19.  Integral  Profile  Matching .  35 

20.  Profile  Matching  with  Integrals .  36 

21.  Six- Point  Profile  Matching .  37 

22.  Profile  Matching  at  Six  Points  Only .  38 

23.  Solid  Angle  Geometry  for  Albedo  and  Earth 

Radiance  Calculation .  41 

24.  Coordinate  Relations  for  Albedo  and  Earth 

Radiance  Calculation .  42 

25.  Albedo-Radiance  Integration  Coordinates .  43 

26.  Area  Elements  for  Cylinder- Cone .  44 


(SCAN  NOTE 


PAGE  IS  BLANK  IN  ORIGINAL  DOCUMENT) 


Earth  gravitational  potential  constant 
Angular  momentum 

Vehicle  rotational  angular  momentum  or  turning  axis 

Orbital  angular  momentum 

Length  of  line  segment 

Mean  anomaly 

Mass  of  earth 

Normal  vector  or  nodal  vector  or  no.  of  line  segments 

Nodal  vector  reference  for  angle  4>a 

Nodal  vector  reference  for  angle  <£v 

Normal  vector  to  vehicle  surface 

Mean  motion 

Vector  projection  of  R  in  plane  perpendicular  to  S, 

Pitch  angle  of  vehicle  relative  to  plane  perpendicular  to  S, 

Distance  from  line  source 

Position  vector  on  earth's  surface 

Position  vector  for  vehicle 

Polar  distance  on  x*  y*  plane 

Earth  radius 

Magnitude  of  Rv 

Segment  spacing  (on  center),  line  source 
Unit  vector  of  solar  position 
Kth  region  of  area  for  integration 
Period  of  motion 
Orbital  period,  Keplerian  orbit 
Time 

Time  from  pericenter 
Velocity  vector 
Magnitude  of  V 

Total  flux  per  unit  length  of  radiating  line  segment 


vm 


AEDC-TDR-63-206 


xa 


xi 

x*  ,y* 
x,y,z 


aN 

P 

y 

6 

OA 

Oe 

Oi 


% 

Oy 

P 

P 

Pe 

a 

T 

<f> 

■<f>A 

0v 

<A 

0 


Coordinate  of  target  for  line  sources 

Coordinate  for  line  segment 

Rectangular  coordinates  in  orbit  plane 

Right-handed,  earth-centered,  inertial,  rectangular 
coordinates 

Azimuth  angle 

a  of  normal  vector 

Nadir  angle 

P  of  normal  vector 

Horizon  angle 

Angle  between  plane  of  motion  and  Sx 
6  from  A  ,  A  plane 
Angle  between  Re  and  S, 

th 

Angle  to  i  line  segment  point,  or  inclination  of 
orbit  plane  to  equatorial  plane 

<?  from  R,  V  plane 

Angle  between  Ry  and  Si 

GMe  ,  constant  of  motion 

Angle  from  normal  for  incident  radiation  at  vehicle 
Roll  angle  of  vehicle 

Distance  between  vehicle  and  point  on  earth's  surface 
Cone  half  angle 

Path  angle  defined  by  Ryt  x  Vt  =  sin  r 
Turning  angle  or  angle  between  Ry,  Re 
Turning  angle  of  A  from  Na 
Turning  angle  of  R  from  N0 

Angle  from  normal  for  reflected  or  emitted  radiation 
at  earth's  surface 

Angle  of  line  of  nodes  from  x  axis  in  equatorial  plane 

Angle  from  line  of  nodes  to  pericenter  in  orbit  plane  or 
solid  angle 


IX 


AEDC-TDR-63-206 


~  Replaces  Q  +  a  when  di  =  0,  since  then  line  of  nodes 
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1.0  INTRODUCTION 


This  report  assembles  the  results  of  studies  performed  by  Scientific 
Computing  Services,  ARO,  Inc.  ,  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  in  the  field  of 
space  flight  thermal  radiation  simulation  to  provide  a  system  to  compute 
real-time  information  about  the  intensity  and  distribution  of  thermal 
radiation  fields  illuminating  a  vehicle  in  near  earth  orbit,  in  transfer 
trajectories,  and  in  orbit  about  any  primary  body.  The  time  functions 
so  obtained  may  eventually  control  the  sequencing  of  parameters  in  the 
Mark  I  aerospace  environment  simulator.  They  will  at  least  provide  an 
independent  means  of  evaluating  performance  of  the  actual  system  when 
it  becomes  operational. 


This  report  is  mainly  limited  to  near  earth  orbits,  but  discussion 
is  sufficiently  generalized  to  indicate  how  transfer  trajectories  and 
alternate  primary  bodies  may  be  incorporated.  We  shall  first  outline 
the  simulation  problem,  then  sketch  a  possible  comprehensive  computer 
program,  and,  finally,  go  into  details  of  the  arithmetic  involved  for  a 
variety  of  options  that  may  be  included.  Appendix  I  is  a  list  of  the 
7074  programs  at  AEDC  concerned  with  Mark  I  illumination  simulation. 


Fig.  1  The  Simulation  Problem 


A  space  vehicle  is  exposed  to  direct  solar  illumination  except  when 
eclipsing  occurs  because  of  shadowing  by  a  primary  body  (Fig.  1).  If 
the  vehicle  is  sufficiently  close  to  the  primary  body,  it  may  receive  sig¬ 
nificant  additional  illumination  of  two  kinds .  The  primary  body  reflects 
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some  fraction  of  the  direct  solar  radiation  incident  on  it  (Albedo)  and 
absorbs  and  re-emits  additional  thermal  radiation  (planetary  radiance). 
We  will  take  the  earth  to  be  the  primary  body  in  what  follows;  the  appli¬ 
cation  to  other  primaries  is  generally  obvious. 

For  a  near  earth  orbit,  the  direct  solar  radiation  intensity  is  taken 
to  be  the  "solar  constant",  invariant  with  orbital  position.  The  direct 
solar  illumination  of  the  vehicle  depends  on  its  attitude  relative  to  the 
sun  and  vanishes  during  eclipse.  The  albedo  radiation  depends  on  effec¬ 
tive  surface  properties  of  the  primary,  assumed  uniform  for  a  particu¬ 
lar  primary,  and  rather  complexly,  on  the  detailed  vehicle-primary-sun 
configuration.  It  vanishes  during  eclipse.  The  earth  radiance  illumina¬ 
tion  depends  on  the  properties  of  the  primary  body,  on  the  vehicle - 
primary  configuration,  and  continues  during  the  eclipse. 

The  paragraph  above  describes  the  illumination  field  in  the  neighbor¬ 
hood  of  an  orbiting  vehicle.  The  thermal  input  to  a  vehicle  depends  upon 
its  shape,  its  orientation  in  the  illumination  field,  and  its  surface  proper¬ 
ties  relative  to  the  spectral  content  of  the  radiation  (Fig.  2). 
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The  illumination  field  becomes  the  thermal  input  to  the  vehicle  at 
its  surface.  Consequently,  the  thermal  input  depends  on  the  particular 
point  on  the  vehicle  (Fig.  3).  Any  surface  element  may  have  distinct 
thermal  properties  and  may  be  exposed  to  all  three,  two,  one,  or  none 
of  the  illumination  components. 


Note:  Surface  properties  do 
not  specifically  apply  to  the 
Agena  vehicle. 


Orbit:  Earth- 1 nertial-Oriented 
or  Sun-Oriented; 
Circular 

Radiation: 

Solar  1400  w/m^ 

Earth  650  w/rrr 

Space  Heat  Sink  -320°F 


b.  Apollo  Type 

Fig.  3  Typical  Shapes  and  Properties 

This  complex  interdependence  is  to  be  simulated  in  the  Mark  I 
chamber,  possibly  synchronized  with  other  simulation  parameters, 
such  as  vacuum  conditions  and  propulsion  loads.  Solar  radiation  will 
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be  simulated  by  the  "solar  bank",  an  extensive  planar  array  of  optical 
elements  to  provide  highly  directional  (collimated)  illumination  of  essen¬ 
tially  true  solar  spectral  composition.  Albedo  and  earth  radiance  are 
to  be  simulated  by  an  open  "bird-cage,  "  essentially  several  linear 
arrays  of  heating  elements  with  roughly  parabolic  cylindrical  reflectors. 
The  solar  bank  is  fixed  in  the  chamber.  The  vehicle  is  mounted  on  a 
horizonal  pitch  axis  so  that  its  pitching  plane  bisects  the  solar  bank. 

The  vehicle  rolls  about  an  axis  that  is  free  in  the  pitch  plane.  The  bird¬ 
cage  array  pitches  with  the  vehicle  but  does  not  roll.  The  whole  range 
of  possible  vehicle-primary-sun  configuration  is  thus  to  be  simulated  by 
pitch  and  roll  of  the  vehicle  and  by  shifting  the  illumination  pattern  of 
the  birdcage. 

Eclipsing  means  simply  extinguishing  the  solar  bank  and  the  albedo 
component  of  the  birdcage  illumination  pattern  but  leaving  the  earth 
radiance  component  turned  on. 

A  further  interdependence  condition  is  imposed  by  the  Mark  I  cham¬ 
ber  itself.  The  system  of  solar  bank,  vehicle  mount,  and  birdcage  is 
contained  in  a  large  vacuum  vessel  whose  dark  refrigerated  inner  walls 
simulate  the  heat  sink  of  space.  The  large  area  of  the  solar  bank  re¬ 
places  some  cold  wall  area,  and  the  solar  bank  is  the  principle  source 
of  the  heat  load  on  the  wall  refrigeration  system.  Consequently,  it  is 
desirable  that  any  element  of  the  solar  bank  shall  be  turned  off  whenever 
its  collimated  beam  will  not  strike  the  vehicle.  The  pattern  of  "on" 
solar  bank  elements  depends  on  vehicle-sun  attitude  and  on  vehicle  shape. 


This  report  will  attempt  to  show  how  the  control  parameters  may 
be  computed  and  inter-related  as  required  by  the  nature  of  any  particu¬ 
lar  test.  The  actual  use  of  the  time-sequenced  parameters  by  the  con¬ 
trol  hardware,  an  instrumentation  problem,  is  beyond  the  scope  of  this 
report. 


2.0  THE  GENERAL  COMPUTER  PROBLEM 

Since  simulation  takes  place  in  real  time,  and  all  control  parameters 
are  to  be  computed  in  advance,  it  is  highly  desirable  that  the  computing 
time  be  a  small  fraction  of  the  corresponding  real  time  for  a  test  of  any 
significant  duration.  This  is  especially  true  if  the  computing  is  done  on 
a  general  purpose  machine  which  has  additional  workloads.  Two  ways 
can  be  exploited  to  make  the  computer-to-real-time  ratio  more  favor¬ 
able.  One  may  plan  to  u.s.e  a  faster  computer  or  to  simplify  the  calcu¬ 
lations,  at  a  reduction  of  realism  in  simulation. 
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At  the  outset,  the  conditions  to  be  simulated  may  be  chosen  so  that 
they  compute  simply.  An  example  is  a  circular  near-earth  orbit  with 
all  parameters  fixed  once  and  for  all,  adequate  for  very  few  revolutions, 
over  a  time  interval  less  than  a  day,  after  which  significant  changes 
must  be  introduced  if  the  simulation  is  not  to  degrade  seriously. 

At  a  more  advanced  level,  for  sophisticated  simulation,  the  com¬ 
puter  time  rises  because  of  the  form  of  equations  which  must  be  evalu¬ 
ated.  Here  it  is  possible  to  introduce  approximating  functions  which 
compute  faster. 

These  considerations  will  be  dealt  with  as  they  arise  in  the  follow¬ 
ing  pages.  In  the  past,  our  work  has  been  based  on  the  idea  that  we 
strive  for  the  highest  realism  until  it  proves  impractical.  We  have 
then  at  least  a  measure  of  the  cost  of  coarser  approximation  in  loss  of 
realism.  The  geometry  of  the  simulation  forces  a  number  of  approxi¬ 
mations;  our  lack  of  certain  detailed  knowledge  and  our  inability  to 
handle  massive  detail  force  further  approximations. 

The  most  feasible  program  appears  to  be  one  consisting  of  two  or 
three  blocks,  as  required  by  the  size  of  the  problem  and  computer 
storage  needed  (see  Table  1). 

In  the  first  block,  all  orbit  calculations  are  performed,  and  from 
them  are  found  the  roll,  pitch,  their  rates,  and  the  existence  of  either 
the  eclipse  condition  or  the  excess  roll  rate  condition  (these  two  cannot 
occur  together  in  principle).  Also,  in  case  of  monotonic  pitch,  revolu¬ 
tions  in  pitch  are  counted. 


TABLE  1 

PROGRAM  BLOCK  OUTLINE 

CALL  BLOCK  1 


Read  input,  initial  conditions,  and  fixed  test  parameters. 

Calculate  orbit  parameters 

Block  1 

Test  eclipse 

Orbit, 

Roll,  pitch,  rates 

Eclipse, 

(Count  revolutions  in  pitch,  if  they  occur) 

Roll, 

Output  tape  A 

Pitch 

Advance  time,  etc.  (introduce  perturbations  ? ) 

Rates 

— Repeat 
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TABLE  1  (Concluded) 

CALL  BLOCK  2 


Read  tape  A,  first  record. 

Block  2 

r— ►  If  eclipsing,  set  up  to  unwind  pitch:  if  multiple 

Housekeeping 

and 

revolutions,  set  earth  radiance  birdcage  pattern. 

Pattern 

If  not  eclipsing,  test  roll  rate. 

Settings 

If  roll  rate  excessive,  smooth  out. 

If  not  eclipsing,  set  solar  bank  pattern  set  "bird¬ 
cage"  pattern,  albedo. 

Write  control  tape  B  (or  Block  3  work). 

Read  next  record. 

—  Repeat. 

CALL  BLOCK  3 


Read  tape  B,  first  record. 

Block  3 

Convert  to  control  media. 

Conversion 

Write  tape  C. 

to  Final  Form 

Read  next  B  record. 

—  Repeat. 

End. 

In  the  second  block,  the  roll  function  is  smoothed  if  it  has  excessive 
rate,  and  the  illumination  array  patterns  are  reset.  During  the  eclipse 
condition,  the  monotonic  pitch  condition  may  indicate  that  pitch  be  backed 
off  one  or  more  revolutions  rapidly  to  unwind  control  cables.  This  con¬ 
dition  occurs  only  when  eclipsing  times  are  maximum.  Part  of  this 
block  may  have  to  be  delayed  to  the  next  block. 

Block  three  is  a  simple  conversion  routine  to  prepare  the  final  con¬ 
trol  media,  either  tape  or  cards,  in  the  format  required  by  the  control 
hardware. 
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Blocks  in  which  computing  predominates  are  built  up  from  program 
elements  for  which  options  are  to  be  selected  for  specific  tests  (see 
Table  2.) 


TABLE  2 

PROGRAM  ELEMENTS 

I.  THE  ORBIT 

A.  Earth  Centered 

1.  circle 

2.  ellipse 

3.  exact 

B.  Perturbations 

4.  diurnal  motion  of  the  sun 

5.  oblate  earth  perturbations 

6.  other  body  perturbations 

C.  Other  Orbits 

7.  orbits  about  other  primaries 

8.  transfer  orbits 

9.  launch- recovery  trajectories 


II. 

SOLAR 

SIMULATION 

A. 

Roll,  Pitch,  Rates 

B. 

Eclipsing 

C. 

Editor  for  Mechanical  Limits 

D. 

Solar  Bank  Pattern  Control 

III. 

SECONDARY  ILLUMINATION 

A. 

Albedo-Earthshine  Field 

B. 

"Birdcage"  Pattern  Control 

3.0  SOLAR  CONTROL  AND  ORBITS 


3.1  ORBIT  DEFINITION 

It  is  intuitively  desirable  to  think  of  an  orbit  as  a  closed  curve,  a 
circle,  or  an  ellipse,  although  this  may  be  far  from  realistic.  We 
define  the  "osculating  ellipse"  as  that  ellipse  which  instantaneously 
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matches  the  actual  orbital  path  and  consider  the  orbit  in  this  manner. 

In  the  absence  of  perturbations,  simple  one  body  central  force  motions, 
the  osculating  ellipse  has  constant  parameters  and  provides  an  ideal 
stationary  description. 

In  the  presence  of  perturbations,  we  may  think  of  the  osculating 
ellipse  as  being  slowly  changed  in  shape,  in  orientation  in  its  plane, 
and  its  plane  slowly  changing  orientation  in  space.  These  changes  are 
mostly  oscillatory;  only  those  caused  by  non -sphericity  of  the  earth 
provide  significant  unidirectional  changes.  We  will  look  at  this  later. 


3.2  ORBIT  CALCULATIONS 

Whatever  mathematical  model  is  used  to  calculate  the  orbital  motion, 
it  must  produce  certain  vector  quantities  or  their  equivalents.  We  re¬ 
quire  the  position  vector  Ry ,  velocity  vector  V  ,  and  perhaps  the  vehicle 

axial  vector  At  ,  which  defines  the  axis  of  the  spin  (roll  axis),  and  its 

— ^ 

time  derivative  A.  (Note:  subscript  1  means  unit  vector.)  (See  Figs.  4 
and  5). 


Fig.  4  Distinction  between  Rotational  Quantities  and  Revolutional  Quantities  in  Flight 
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Lv 


In  any  plane  circle  there  is 
at  least  one  diameter  (N,  l\T) 
perpendicular  to  any  given 
direction  Sj. 

The  positive  sense  of_ vector  oN 
is  taken  as  that  of  (Sj  x  lq)p 

The  angle  I  defining  Rj  is 
measured  from  N. 

Fig.  5  Calculation  of  Node  Vectors 

For  any  orbital  motion,  the  first  immediately  useful  information 
obtainable  is  the  eclipsing  interval.  Eclipsing  depends  only  on  the  posi¬ 
tion,  not  the  attitude,  of  the  orbiting  body  relative  to  the  earth's  shadow. 
The  earth  is  assumed  to  be  a  sphere,  devoid  of  atmosphere,  casting  a 
right  circular  cylindrical  shadow,  geometrically  defined  (Fig.  11,  page  19). 

— > 

If  Ry  is  the  position  vector  of  the  orbiting  vehicle,  Re  ,  the  earth 
radius,  and  is  a  unit  vector  from  earth  to  sun,  we  find  the  projec¬ 
tion  P  of  Ry  on  a  plane  perpendicular  to  Sx  ,  and  test  the  magnitude  of 
the  projection  against  Re  .  If  |  P  |  <  Re  and  if  the  vehicle  is  on  the  aphelion 

limb  of  the  orbit,  the  vehicle  is  in  eclipse.  The  last  condition  is  more 

— >  — > 

simply  stated  as  the  requirement  that  Ry  •  S,  be  negative.  The  depen¬ 
dence  of  eclipse  interval  on  orbit  parameters  is  detailed  in  Ref.  1. 

The  eclipse  indicator  may  be  used  later  to  bypass  all  albedo  calcu¬ 
lations,  to  turn  off  the  solar  bank  entirely,  and  to  relieve  the  constraint 
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on  number  of  revolutions  in  pitch  angle.  This  last  situation  occurs  for 
orbits  where  the  vehicle  roll  axis  A1  passes  through  the  sun,  for  when 
this  happens  the  pitch  is  monotonic  rather  than  oscillatory,  as  will  be 
seen  below. 


The  vector  product  Ry  x  V  defines  the  orbital  angular  momentum 
vector  (per  unit  mass)  L  and  hence  defines  the  instantaneous  plane  of 

the  orbit.  The  vector  product  JA,  x  A  )1  defines  the  axis  about  which 
the  vehicle  is  turning  its  axis  A  at  the  instant,  and  this  is  the  vehicle 
angular  momentum  vector  L A  (if  there  is  no  spin  about  A  ). 

— 4  — > 

For  some  modes  of  orbital  flight,  the  two  vectors,  L  and  La,  coin¬ 
cide.  Hence  the  turning  axis  (in  the  direction  of  La)  may  be  obtained 
directly  from  the  orbital  vectors  R  and  V  .  For  other  modes  of  flight. 
La  may  not  be  determinate  from  the  orbit  but  must  be  specially  given. 

—4 

Now  if  S,  is  a  unit  vector  pointing  to  the  sun,  the  angle  0 s  between 
orbit  plane  and  sun  direction  is  obtained  from 

-4 

St  •  Lj  m  —  sin  (9g 

and  a  nodal  vector  from  which  orbital  angles  are  measured  is 

N0  =  (S,  X  Ll)i  (Fig.  5) 

The  angle  <£v  between  Ry  and  N0  is  then  obtained  from 

->  -4 

COS  =  Ryt  •  No 


Similarly,  at  any  instant  the  vehicle  axis  vector  A  is  turning  about 
the  vector  La  in  a  plane  inclined  at  angle  9 a  defined  by  Si  •  La  =  -  sin  9 a 
to  the  sun  direction  S, .  A,  has  turned  through  angle  4>a  from  a  nodal 

-4  —4—4 

direction  Na  =  (S,  x  La )t  such  that 

-4  —4 

cos  (f) a  =  At  •  Na  • 

“4  -4 

As  defined,  the  nodal  vectors  N0  and  Na  lie  in  their  respective  turning 
planes  and  jire  perpendicular  to  the  sun  vector  as  well  as  to  their 
respective  L  and  La-  If  the  angular  momentum  L  and  La  always  coin- 

,  —4  -4  —4  —4 

cide  in  direction,  then  Ry,  V,  A  and  A  are  coplanar,  the  node  vectors 
N0  and  Na  coincide  and  A  and  Rv  are  related  by 

0 A  =  0V  +  T 


where 


tan  r 


R  V  X  A 

~4  -4 

Ry  .  a 
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The  distinction  between  the  pairs  of  vectors  Ry  ;  V  and  A,  A  is  as 
follows:  Ry  and  V  are  orbit  descriptors,  serving  to  determine  vehicle 
location  in  orbit,  essential  in  establishing  the  eclipse  condition,  and  for 
locating  the  vehicle  position  relative  to  the  Albedo  source,  the  sunlit 

— > 

hemisphere  of  the  earth.  The  vectors  A  and  A  are  vehicle  attitude 
descriptors,  governing  its  exposure  to  direct  solar  radiation  and  its 
attitude  toward  the  earth  for  both  Albedo  and  earthshine  radiation.  For 
some  orientations  of  the  vehicle  in  orbit,  these  vector  sets  are  func¬ 
tionally  coupled,  and  only  one  needs  to  be  developed. 


3.3  COORDINATE  SYSTEM 

Although  we  use  vector  notation  for  freedom  from  particular  co¬ 
ordinate  systems,  numerical  computations  do  require  vector  components 
which  are  dependent  on  the  coordinate  system  selected. 

In  the  following,  the  fundamental  coordinates  are  geocentric,  non¬ 
rotating,  orthogonal  axes  defined  as  follows.  The  z-axis  points  north, 
parallel  to  the  earth's  axis.  The  x-axis  points  to  the  vernal  equinox  of 
some  particular  date,  *  and  the  y-axis  is  so  located  as  to  form  a  right 
handed  system  (Fig.  6).  As  required  for  clarity  or  efficiency,  other  co¬ 
ordinate  systems  will  be  defined  and  introduced. 


N 

Z 


Vernal  Equinox 
(of  standard  date) 

Fig.  6  Earth-Centered  Inertial  Coordinate  System 


*The  date  of  reference  for  particular  tables;  for  example,  "the 
mean  equinox  of  1950.  0.  " 
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We  use  the  geocentric  inertial  system  so  that  astronomical  data 
may  be  used,  directly  from  "Ephemeris"  tables  (Ref.  2).  For  any 
degree  of  realism  in  computing  attitude  of  a  vehicle  to  the  sun,  we  re¬ 
quire  coordinates  of  the  sun  position  relative  to  the  earth.  For  high 
realism  other  body  coordinates  will  be  needed.  These  are  available  on 
"Themis"  magnetic  tapes  (Ref.  3). 

Within  the  Mark  I  chamber  the  coordinates  for  the  motion  are 
pitch  p  and  roll  p ,  angles  defined  by  the  mechanical  construction  as 
previously  described.  Depending  upon  the  constraints  of  the  vehicle 
motion,  the  sun  angle  6  and  vehicle  turning  angle  <f>  will  be  defined  in 
slightly  differing  terms,  and  from  these  angles,  roll  and  pitch  are 
calculated. 

3.4  ORIENTATION 

By  vehicle  orientation,  we  mean  certain  fixed  constraints  of  rela¬ 
tions  between  the  vehicle  axes  and  directions  in  space  or  to  other  bodies. 
There  are  five  commonly  conceived  orientations  (Fig.  7): 

1.  A  vehicle  which  has  no  motion  about  its  center  of  mass 
is  space  oriented.  It  may  have  only  one  axis  fixed  in 
direction,  typically  the  spin  axis  for  a  spin  stabilized 
missile. 

2.  A  vehicle  may  have  a  fixed  orientation  relative  to  the 
earth,  for  example,  a  photo-reconnaisance  vehicle,  and 
is  said  to  be  earth  oriented. 

3.  A  minimum  of  solar  radiation  input  to  an  orbiting  vehicle 
is  obtained  when  minimum  cross  section  area  is  exposed 
to  the  sun,  a  solar  orientation. 

4.  The  roll  axis  of  the  vehicle  may  be  maintained  coincident 
with  the  velocity  vector,  tangent  to  the  orbit,  thus  called 
tangent  oriented.  This  is  the  most  common  concept  of 
orbital  flight  orientation. 

5.  If  the  orbital  angular  velocity  vector  and  the  vehicle 
angular  velocity  vector  coincide  and  are  equal  (i.  e. ,  the 
rotation  period  equals  the  revolution  period  and  the 
motions  are  coaxial),  this  may  be  called  an  angular 
momentum  orientation.  This  will  be  also  tangent  orienta¬ 
tion  for  a  circular  orbit  and  will  be  earth  oriented. 
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Space 

Oriented 

Invariant  Direction 
in  Space 


Earth 

Oriented 


Quarter  Period 
Segments 


Constant 

Rotational 

Angular 

Momentum 


Fig.  7  Orientation  Types 
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From  these  examples  we  see  that  orientation  can  be  partly  or  com¬ 
pletely  specified  by  stating  the  fixed  relation  between  one  or  more  axes 
of  the  vehicle  and  a  fixed  space  direction,  a  direction  to  another  body, 
or  relative  to  definite  directions,  in  one,  two,  or  three  dimensions.  A 
vehicle  employing  thrust  for  attitude  control  may  have  continually  vary¬ 
ing  orientation. 

If  we  are  to  obtain  analytic  descriptions  of  orbital  flight,  it  is 
desirable  to  assume  a  definite  orientation  to  hold  throughout  the  time 
interval  of  interest.  In  the  following,  we  shall  assume  tangent  orienta¬ 
tion  without  spin.  We  will  remain  aware  of  the  fact  that  a  real  flight 
with  attitude  control  may  be  piecewise  analytic  but  not  describable  in 
advance,  unless  scheduled  accurately.  In  that  case,  the  ordinary  laws 
of  mechanics  can  be  applied  to  determine  the  complete  history  of  orien¬ 
tations  in  advance. 

Tangent  orientation  means  that  the  vehicle  axis  lies  in  the  plane  of 
the  orbit.  For  a  circular  orbit,  the  axis  vector  A  will  be  perpendicular 
to  the  position  vector  Ry.  Hence  </>A will  be  90  deg  ahead  of  <£v  •  For 
elliptical  orbits,  the  axis  vector  will  be  perpendicular  to  the  position 
vector  only  at  perigee  and  apogee.  From  perigee  to  apogee  the  vehicle 
is  gaining  altitude,  the  axis  is  above  the  horizonal;  during  the  passage 
apogee  to  perigee,  it  is  below  the  horizonal.  For  ordinary  Keplerian 
orbital  motion,  this  angle  of  elevation^or  depression  is  easily  calcu¬ 
lated,  so  that  the  vehicle  axis  vector  A  may  be  found  from  the  orbit 
position  vector  Ry. 

Regardless  of  orientation,  it  is  always  necessary  to  compute  the 
progress  of  the  vehicle  along  the  orbit,  if  only  to  test  for  and  simulate 
eclipsing.  The  orbit  position  also  enters  Albedo  calculations. 

3.5  ROLL  AND  PITCH 

We  shall  calculate  the  Mark  I  control  parameters  roll  and  pitch  from 
the  angles  and  0A.  At  the  outset  we  see  that  the  vector  A  corresponds, 
not  to  a  particular  axis  in  the  vehicle,  but  to  an  axis  determined  by  the 
vehicle  mounting  in  the  Mark  I  cell.  The  vector  A  is  the  cell-defined  roll 
axis,  and  obviously  the  reference  direction  for  roll  can  be  anywhere  in 
the  vehicle,  as  long  as  it  is  perpendicular  to  A  .  We  assume  the  conven¬ 
tion  that  the  zero  pitch  and  roll  position  is  the  original  mounting  configu¬ 
ration  in  the  cell.  The  initial  pitch  and  roll  for  the  motion  depend  on  the 
earth-sun-vehicle  configuration  at  injection  into  the  orbit. 

A  vehicle  on  which  no  torques  act  has  constant  rotational  angular 
momentum  La  (Fig.  8).  If  this  vector  is  aligned  with  the  sun  direction. 
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Fig.  8  Roll  and  Pitch  Geometry 


La,  =  -  S,  ,  the  plane  of  the  motion  is  perpendicular  to  the  sun  direction. 
Roll  and  pitch  are  constant  in  this  case,  and  it  is  convenient  to  take  this 
configuration  as  defining  the  zeros  of  roll  and  pitch.  This  corresponds 
to  a  circumpolar  orbit,  tangent  oriented  vehicle,  crossing  north  over 
the  equator  at  sunrise.  By  the  definition  of  solar  inclination  angle  <9 a  , 
this  case  means  0a  =  n/2  .  If  we  let  6\  =  0  ,  the  plane  of  the  motion  con¬ 
tains  the  sun  vector  S,.  Any  time  the  vector  A,  coincides  with  +S,,  the 
roll  p  is  undefined.  At  such  critical  points  we  have  an  operating  option: 
we  may  allow  pitch  p  to  increase  monotonically;  or  we  may  reverse  the 
pitch  direction  and  change  the  roll  by  180°.  In  the  first  option  roll  is  un¬ 
changed,  and  pitch  can  go  through  several  complete  revolutions  before 
mechanical  restrictions  force  us  to  "unwind".  In  the  second  option,  roll 
can  be  flipped  only  as  fast  as  the  roll  drive  permits  (about  1  rpm),  but 
pitch  is  then  limited  to  the  range  -  n/2  ,  0  ,  n/2  ,  and  there  is  no  unwind¬ 
ing  problem.  All  this  applies  only  for  0a  =  0  .  If  the  sense  of  positive 
pitch  is  toward  the  sun,  with  zeroes  as  described  above,  we  have 
pitch  p  and  roll  p  given  by  spherical  trigonometry  in  the  forms 

sin  p  =  cos  0a  s'n  S^A 

sin  p  =  cos  6a  cos  ^/cos  p. 

From  these  expressions  we  see  that  sin  p  is  limited  in  amplitude  to 
cos  O a  an<j  is  otherwise  simply  proportional  to  sin  <f>A  as  <f>A  goes 
through  a  complete  revolution,  p  oscillates  smoothly  between  ±  ( n/2  -  B a)  • 
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If  <9a  =  n/2  ,  the  oscillation  has  zero  amplitude,  i.  e.,  p  =  0  .  Finally  if 
dA  =  0  ,  we  have  simply  p  =  0a  or,  optionally 

For  0A  =  0,  ±  n/2  ,  p  =  0  dp  =  d  0 

0  A  =7 7  +  77/2  p=77-0  dp=-d0 

Examining  the  expression  for  p  in  the  same  way,  we  see  p  is  oscilla¬ 
tory  in  0a  with  amplitude  determined  by  dA .  We  note  that  dA  =  n/2  gives 
cos  p  cos  p  =  1  ,  requiring  p  =  0  andp  =  0  ,  and  <9a  =  0  gives  sin  p  =  1  , 

P  =  77/2  .  At  intermediate  0A,  p  =0  gives  sin  p  =  cos  0A  cos  0a.  When 

both  are  oscillatory,  p  is  n/2  out  of  phase  with  p  .  When  0A  is  zero,  and 
we  choose  the  limited  pitch  option  above,  we  have  for  p  these  conditions: 

when  0  <  0a  <  77 ,  p  =  n/2 

77  <  0a  <  2  77  ,  p  =  -  77/2 

These  relations  are  also  implied  in 

cos  p  cos  p  =  sin  0a 

for  when  0A  vanished,  since  p  is  simply  0A  or  77  -  0A  ..except  for  0A  =  ±  n/2  , 

we  must  have  p  identically  ±  n/2  .  And  when  dA  =  n/2  ,  we  must  have  both  P 

and  p  =  0  or  77 . 

We  now  examine  the  roll  and  pitch  rates  p  and  p  .  By  differentiating 
and  substituting  identities  we  find 

P  =  sin  p  0A 
p  =  —  tan  p  cos  p  <pA 

where  0A  is  the  instantaneous  angular  velocity  of  rotation  of  vector  A 

-> 

about  its  axis  La  .  This  form  ignores  any  contribution  caused  by  pre¬ 
cession  which  would  give  terms  in  dA  .  Such  terms  are  normally  much 
smaller  than  those  in  0a. 

We  see  that  p  has  a  singularity  if  p  =  ±  77/2  .  From  earlier  argu¬ 
ments,  this  can  only  occur  when  6 a  =  0  ,  but  for  dA  =  0  ,  p  =  ±  n/2  ,  cos  p, 
vanishes  and  p  is  indeterminate.  For  <?a  small  p  takes  on  large  values 
at  p  =  ±  -y  (Fig.  9).  The  Mark  I  mechanical  limit  of  1  rpm  means  that 
we  smooth  across  the  steep  portion  of  the  graph  of  roll  versus  0A  »  with 
a  line  or  curve  having  slope  not  exceeding  1  rpm.  Except  in  extreme 
instances,  less  time  is  involved  in  this  smoothed  portion  if  we  delay 
smoothing  until  required,  rather  than  try  to  anticipate  it  (Fig.  10),  and 
the  computer  programming  is  considerably  simplified. 
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This  condition,  roll  rate  exceeding  the  mechanical  limit,  must  be 
identified  when  computing  roll  and  pitch.  The  infinite  values  do  not 
actually  occur  since  they  are  absorbed  by  the  pitch  becoming  monotonic 
rather  than  oscillatory,  or  optionally,  by  reversing  the  sign  of  p  as 
stated  earlier  and  flipping  p  through  n  ,  which,  of  course,  is  also  subject 
to  the  1-rpm  limit.  At  worst,  only  1/2  min  suffices  to  switch  over. 

This  can  occur  twice  per  revolution  in  </>A  »  and  its  importance  is  thus 
dependent  on  .  For  tangent -oriented  orbital  flight  simulation,  we 
have  the  excess  rate  condition  occupying  one  minute  or  less  out  of  the 
orbit  period  of  (nominally)  90  min  or  more. 

The  smoothing  of  the  roll  function  to  avoid  excess  roll  rate  will  be 
done  by  an  editing  program  block.  As  the  flight  advances,  roll,  pitch, 
and  rates  are  computed,  roll  rate  is  tested  and,  if  excessive,  an  indi¬ 
cator  is  turned  on.  All  this  occurs  in  the  first  program  block.  The 
editing  program  block  will  scan  for  the  excess  roll  rate  indicator  and, 
when  it  is  found,  will  compute  roll  values  from  the  smoothing  line  until 
the  rate  is  again  within  the  allowed  range. 

In  the  next  section,  we  shall  follow  the  calculation  for  the  simplest 
kind  of  orbital  flight  model. 


3.6  ORBITS 

3.6.1  Circular  Orbit,  No  Perturbations 


In  uniform  rotational  motion,  we  have  the  progress  around  the 
orbit,  through  central  angle  <£y  » 

J.  jl  .  2  77-  (  t  -  t„  ) 

9V  =  9Vo  +  - - 


where  4>yo  is  the  initial  angle,  t0  is  the  initial  time,  Ty  is  the  orbit 
period. 


Ty  =  2  7T 


with  rv  the  orbit  radius,  G  the  universal  gravitation  constant,  and  Me 
the  mass  of  the  central  body. 

We  assume  tangent  orientation;  hence  the  path  angle  r  =  n/2  ,  and 
the  roll  axis  vector  A  turns  through  angle  <f>A  , 

</>  A  =  </>V  +  n/2 

The  node  vectors  coincide,  NA  =  Ny  ,  as  do  the  angular  momenta, 
Ly1  =  LAl  and  the  planes  of  the  motions  so  that  0y  =  dA .  We  assume 
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that  injection  into  orbit  (time  to)  occurs  at  sunrise  at  a  latitude  equal 
the  sun  declination  of  date,  and  hence,  at  the  nodal  point,  and  at  an 
altitude  h  such  that  h  =  rv  -  re  (  re  is  earth  radius).  We  define  the 
horizon  angle  y  by 

sin  y  =  re/rv 

an  important  quantity  in  secondary  illumination.  Then  eclipsing  occurs 
between  4>y  points  (Fig.  11): 

3  rr/2  +  A^>  if  0  <  ]  0y  |  <  y 


or 


77/2  +  A  if  7r-y<|0y|<77 


and 


A  0 


,in  y  —  sin  Oy 
cos  6tj 


This  is  equivalent  to  the  earlier  test  expressions 

Ry2  -  (  Ry  •  Sl)  <  Re2  and  Ry  •  S,  <  0 

required  for  eclipsing  (Fig.  12). 


-Sun 
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Eclipsing 

b.  View  from  Sun 
Fig.  12  Eclipse  Test 

In  the  circular  orbit,  the  fraction  of  time  in  eclipse  is  plainly  A $/n  , 
which  attains  a  maximum  of  yht  when  the  orbit  plane  contains  the  sun 
direction  ,  Oa  =  0  . 

The  pitch  and  roll  and  their  rates  are  calculated  by  the  equations 
of  the  the  previous  section. 

We  may  evaluate  the  effect  of  excessive  roll  rate  by  the  inequality 
|  p  |  =  |  tan  p  cos  p  <fc  |v  £  1  rpm  =  2  77/ 60  rdn  /  sec 

and  the  definitions 

<£y  =  2rr/Tv 


which  reduces  to  the  problem  of  evaluating  the  relations  between 
<f>  »  0  ,  y  in 

cos2  6  sin2  cf>  sin2  0  sin3  y  ^  re3 

(  1  —  cos2  6  sin2  cf>  )2  (60)2  G  Me 

where  the  right  side  is  a  constant.  Since  for  excess  roll  to  occur  we, 
must  be  near  <f>  =  0  (or  n)  ,  6  =  0  (or  it)  and  the  orbital  velocity  (hence  <t>) 
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is  greatest  for  low  altitude  or  y  near  77/2,  we  can  approximate  the  rela¬ 
tion  for  small  6  and  4>  by 


**  * 


J 


G  Me 


Given  small  6  for  the  particular  orbit  (circular,  tangent  oriented),  when¬ 
ever  (f>Y  is  within  the  intervals 


A  cf>  «  + 


on  either  side  of  </4y  =  0  or  n ,  we  have  the  excess  roll  rate  condition. 


3.6.2  Keplerian  Orbit,  No  Perturbation 

Planetary  motion  in  a  pure  central  force  field  traces  out  an  ellipse 
with  the  force  center  at  one  focus.  Given  the  initial  position  and  velocity 
vectors  R0  ,  V0  of  an  earth  satellite,  the  entire  motion  is  worked  out  as 
follows:  The  angular  momentum  is  conserved;  R„  x  V0  =  L0  =  R  x  V  is 
a  constant  vector.  The  energy  is  conserved,  and  if  a0  is  the  semi¬ 
major  axis  and  fi  =  G  Me  ,  we  obtain  a0  from  the  energy  equation: 

,  -»  ,2 

1  =  2  _  |  V0  I 

ao  I  R  I  n 

The  plane  of  the  orbit  contains  R0  and  V0 ,  perpendicular  to  L0  . 

From  R0  and  V0  ,  we  are  to  find  the  inclination  of  6\  of  the  orbit  plane 
to  the  (x,  y)  plane,  the  orientation  <a0  of  the  semi-major  axis  (specifi¬ 
cally  the  perigee  direction)  in  the  orbit  plane,  and  the  angle  120  between 
the  line  of  nodes  (ascending  node)  in  the  (x,  y)  plane  and  the  x-axis,  as 
well  as  the  angle  0O  from  perigee  to  the  present  position  R0  (Figs.  13 
and  14).  The  order  of  computing  is  as  follows: 

Inclination  - sin"1  [L0Z/|L0|]  (0  <  dio  <  n) 

Longitude  of  ascending  node 

(If  sin  6{0  =  0  ,  set  =  0  )  tan  fi0  =  j  Lox  !/{  -  Loy  } 


Eccentricity  eQ  =  1  -  I  L°  1* 

v  M  aO 

Eccentric  anomaly  E0  =  y-  -  sin 


(0  <  Q0  <  2n) 
(0  <  e0  <  1) 

I  Ro 
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The  sign  of  sin  (E0)  is  the  same  as  the  sign  of  R0  •  V0  •  If  R0  •  Vn  =  0  , 
then  E0  =  0  ;  if  I  R0  I  <  a0  ,  Ec  -  n  if  I  Ro  |  >  a0  ;  if  |  R0  |  =  aQ  and  R0  •  VQ  =  0  , 
we  may  setE0  =  0,  eQ  =  0  and  circular  orbit  calculations  result.  Angle 
from  perigee  (true  anomaly) 


<f>o 


tan 


(I  -  e02/2  sin  E0 
i  cos  En  —  e„  i 


or 


<£0  =  0  if  eQ  =  0. 


Argument  of  latitude 
if  <9io  =  0 

fi0  =  tan 


1  it o z  /sin  0} 0  I 


if  0in  =  0 


fi0  =  tan 


Rox  cos  +  Roy  sin 


1  oy 


Argument  of  perigee 

&>0  =  0  —  4>c 

Note:  The  symbol  tan 


u 


indicates  that  signs  of 


both  terms  are  considered  to  place  angle  in 
proper  quadrant. 


We  now  have  rotation  angles  Ro  >  ^io  >  <^o  of  an  (inverse)  Euler  rota¬ 
tion  from  the  plane  of  the  ellipse  (x*,  y* )  to  the  (x,  y,  z)  coordinate  sys¬ 
tem.  In(x*,  y*),  we  have  x0*  =  r0  cos  9 So,  y0*  =  rosin0o,  r0  =  I  R0  I  . 

-)  — > 

In  other  words,  corresponding  to  the  initial  R0,  V0  in  U,  y,  z)  and 
the  definite  orbit  they  prescribe,  there  exists  one  ellipse  of  the  same 
eccentricity  and  semi-major  axis  in  (x*,  y*)  where  the  perigee  is  aligned 
with  the  x*  -  axis  .  JThe  angle  from  <f>0  is  preserved  in  both  spaces,  as  is 
angle  t o  between  R0,  VQ. 


We  will  follow  the  motion  in  the  (x*  y*)  plane,  and,  through  the 

— >  —> 

Euler  rotation,  recover  the  corresponding  R,  V  in  the  (x,  y,  z)  space. 
Then  at  anytime,  having' R,  V,  S,in  U,  y,  z),  we  may,  assuming 
tangent  orientation,  compute  0A,  &A  from  equations  given  earlier; 
hence,  roll  and  pitch.  For  initial  roll  and  pitch  rates,  we  need  <£Ao  : 

•  • 

^Ao  “  <£o  A  To 


T 


o 


1  Lp  1  /  (l 

I  R0  I'  \  I  R0  M  V0  I 
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hence 

where 

Mean  anomaly; 


$ Ao  — 


1  L0  |  /x 


->  ,3  ,  ,2 

R„  Vj 


fX  sin  r0 

I  R„  I  I  I 


^0  =  --  -  sin-1  (R0l  •  v0l) 
Mo  =  E0  -  e o  sin  E0 


Time  since  perigee 

tp o  ~  Mo/no 

Epoch  of  perigee  passage 

tp  =  t0  -  tp o 


n  o 


(Keplers  equation) 

—  (mean  motion) 

lo 

( t0  is  initial  time  of  R0,  V0. ) 


To  advance  the  orbit  by  a  time  increment  T,  we  advance  the  time 
since  perigee 

t  =  tpo  +  At  =  t0  -  tp  +  At 

then  M(t)  =  n0(tpo  +  At)  =  n0(t) 

and  E(t)  -  e0  sinE(t)  =  M(t) 

is  solved  by  iteration  for  E(tk 


<£(t)  =  tan 


{  (1  -  e02)1/2  sin  E(t)  i 


r(t) 
x*  (t) 
y*  (t) 


{  cos  E  ( t )  —  e0 
a0  (1  -  eQ  cos  cf>{  t)  ) 
r(t)  cos  0(t) 
r(t)  sin  0(t) 


r(t)  =  J/2 


r(t)  =  sin 


1 


(t) 


r(t)  V(t) 


for  sin  <p  ( t )  >  0 


■(t)  =  n  —  sin  1 


(  r  ( t  TV  Itl)  for  Sin<^(t)  <  0 


t{  t) 
v*(t) 
Vy(t) 

^A(t) 


TT 

2 


for  sin  <f>  (t )  =  0 


v(t)  cos  ( c£(t )  +  r(t)  ) 

v(t)  sin  (</>(t)  +  r(t)  ) 

/ i  sin2  r(t) 
r(t)  |  L0  | 
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We  now  have  new  vectors  R(t),  V  (t)  in  (x*,  v*)  and  can  apply  the  in- 
verse  Euler  rotation  to  each  using  fl0,  tdo,  «0  to  obtain  (R(t),  V(.t))  com¬ 
ponents  in  (x,  y,  z).  From  these,  as  before,  we  develop  roll  and  pitch 
values:  using  cf>A  ,  we  then  obtain  their  rates. 

3.6.3  Keplerian  Orbit  with  Perturbations 


All  the  work  of  the  previous  section  holds  with  the  exception  of  the 
Euler  rotation  from  the  orbit  plane  to  the  tx,  y,  z)  space.  The  inclina¬ 
tion  Oi  of  the  orbit  plane  remains  unchanged,  but  the  perigee  advances 
or  retrogresses,  and  the  plane  of  the  orbit  precesses.  The  earth's 
oblateness  is  the  cause,  through  modification  of  the  earth's  potential 
(Fig.  15). 


Advance  of  Apsidal  Line 
in  Orbit  Plane 

(Change  of  u) 


of  Nodes) 
(Change  of  Q) 


Combined  Effects 


Fig.  15  Secular  Perturbations  Caused  by  Earth  Oblateness 


At  any  time,  withR,  V  found  in  the  (x*,  y*)  space,  having  B  =  JRe2 , 
(J  is  a  gravitational  potential  constant),  we  calculate 

_  _  B  (  cos  0;  )  (  t  —  t0  )  n0  B  (  2  -  2.5  sin2  )  (  t  -  t0  )  n„ 

12  —  120  ~  I  2  ^  2  .  2  \  2 

a02  (1  -  e2)2  a02  (  1  -  e2) 

nQ  is  the  mean  motion  calculated  at  time  tQ. 

<PQ,  fl0,  and  co0  are  found  immediately  from  R,  V  as  given  before.  We 
assume  that  e,  a,  and  di  are  constant,  although  they  are  actually  oscil¬ 
latory. 
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The  Euler  angles  0,  a  are  updated,  then  the  rotation  is  invoked  to 
generate  the  R,  V  components  in  (x,  y,  z)  space,  and  we  proceed  as 
before.  This  process  corresponds  in  part  to  the  Hansen  method  of  com¬ 
puting  satellite  orbits  (Ref.  4). 

We  believe  that  this  approach  to  a  simulation  orbit  is  novel.  It  is 
important  since  the  secular  variations  may  amount  to  as  much  as  5  deg 
per  day  for  near  earth  satellites  (Ref.  5). 

3.6.4  “Exact”  Orbits 

So  called  exact  orbits  are  obtained  as  solutions  to  the  differential 
equations  of  motion, 

miv  =  SF  =  -  l  -GrMA  Rvi  +  SFj 

v  Hyj  )  j 

in  which  the  sum  in  i  is  taken  over  all  major  bodies  concerned  with  the 
motion  as  if  each  were  a  simple  central  force  mass  (assumed  homo¬ 
geneous,  spherical,  or  point),  and  the  sum  in  j  is  over  all  other  forces, 
whatever  their  origin,  which  must  be  known  as,  generally,  functions  of 
position  and  time. 

A  numerical  solution  is  obtained  by  marching  ahead,  keeping  all 
force  terms  up  to  date,  from  the  initial  coordinates  and  initial  time. 
Several  methods  are  available  for  this,  all  requiring  fine  increments 
and/or  multiple  precision  arithmetic,  all  slow  to  compute.  We  have 
found  that  high  precision  orbits  are  costly  and  may  easily  require  "real 
time"  to  compute;  hence,  they  are  almost  useless  for  long  time  schedul¬ 
ing  of  simulation  runs.  (These  tests  were  programmed  for  an  IBM  7070). 
They  are  however  necessary  for  evaluating  any  approximate  method,  such 
as  that  of  the  previous  section. 

These  are  two  historical  numerical  methods  of  solution  of  the  equa¬ 
tions  of  motion.  The  method  of  Cowell  is  quite  straightforward,  in  con¬ 
junction  with  a  Runge-Kutta  numerical  technique.  It  involves  small  dif¬ 
ferences  in  large  numbers  and  requires  multiple  precision  computer 
techniques. 

The  method  of  Enke  is  a  perturbation  method  in  which  the  initial 
conditions  establish  an  osculating  ellipse.  Small  departures  from  the 
osculating  ellipse  are  calculated  until  they  grow  excessive;  then  a  new 
osculating  ellipse  is  calculated.  The  method  yields  good  results  without 
resorting  to  multiple  precision. 

For  each  of  these  methods,  the  inclusion  of  other  planetary  bodies 
requires  reference  to  tables  of  positions,  the  ephemerides,  which  are 
available  on  copies  of  the  "Themis"  tapes  (Ref.  3). 
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Modern  methods  involve  substitution  of  some  faster  techniques  than 
the  Runge-Kutta  procedure  in  solving  the  equations.  These  methods  are 
not  self  starting  but  may  depend  on  a  Runge-Kutta  beginning,  after  which 
they  require  fewer  evaluations  of  the  force  terms  in  the  differential 
equations  (Ref.  4). 

If  no  velocity  dependent  forces  are  involved,  the  principle  secular 
perturbations  are  those  previously  described  caused  by  earth  oblateness, 
all  other  perturbations  being  primarily  oscillatory.  We  believe  that  the 
method  proposed,  a  single  Keplerian  osculating  ellipse  with  secular 
perturbations  entering  an  Euler  rotation,  is  adequate  for  realistic  simu¬ 
lation  and  is  fast  to  compute. 

3.6.5  Arbitrary  Flight  Paths 

It  may  be  required  at  times  to  simulate  non-periodic  motions  of 
vehicles;  for  example,  launch,  recovery  paths,  and  transfer  orbits. 

The  most  general  method  for  this  is  the  direct  numerical  solution  of  the 
differential  equations  of  the  motion,  for  which  a  complete  schedule  of 
thrust  magnitude  and  direction  will  be  required.  For  approximation 
purposes,  a  design  trajectory  may  be  adequate  with  the  position  vector 
given  as  a  function  of  time.  From  this,  the  velocity  vector  is  comput¬ 
able.  An  assumption  of  vehicle  orientation  relative  to  the  trajectory 
allows  complete  motion  description,  and  when  the  initial  date  (accurate 
to  better  than  the  simulation  increment,  say  two  minutes)  is  given,  roll 
and  pitch  can  then  be  calculated  by  previous  methods. 

If  we  are  concerned  only  with  solar  radiation,  only  the  coordinates 
of  the  vehicle  axis  (vector  A  of  earlier  sections)  need  be  specified.  The 
actual  coordinates  of  position  are  needed  if  we  are  interested  in  phe¬ 
nomena  caused  by  another  body,  say  earth,  such  as  eclipsing  or  Albedo 
illumination. 

3.7  SOLAR  BANK  CONTROL 

Only  those  solar  bank  modules  are  to  be  on  which  contribute  to  the 
illumination  of  the  vehicle.  All  others  are  to  be  off  (see  section  1.  0). 
Since  the  vehicle  pitches  relative  to  the  sun  bank,  we  may  consider  the 
projection  of  the  vehicle  outline  onto  the  solar  bank.  Then  any  solar 
module  entirely  outside  the  outline  is  to  be  off,  any  entirely  within  shall 
be  on.  What  about  the  intermediate  cases,  partly  within,  partly  without? 

We  provide  three  options  -  a  module  is  inside  the  outline  (on) 

1.  if  any  part  of  it  is  inside, 

2.  if  its  center  is  inside, 

3.  if  it  is  entirely  inside. 


27 


AEDC-TDR -63-206 


In  addition,  if  the  vehicle  is  in  eclipse  behind  the  earth,  all  modules 
are  off.  There  will  be  no  intermediate  cases  for  eclipse.  An  option  is 
provided  to  either  skip  the  outline  calculations  or  to  continue  them  to 
show  roughly  the  vehicle  outline. 

The  problem  is  general  as  long  as  the  vehicle  consists  of  the  follow¬ 
ing  elements: 

1.  a  spherical  nose, 

2.  a  cone  frustum  tangent  to  the  spherical  nose,  and 

3.  a  right  circular  cylindrical  element  of  same  size  as  the  larger 
base  of  the  frustum. 

Options  allow  the  length  of  any  of  these  to  be  zero.  Hence  most 
combinations  are  available  with  the  restrictions  that  the  order  be  pre¬ 
served  and  the  sphere  diameter  be  no  larger  than  the  cylinder  diameter. 
The  body  must  be  everywhere  convex. 

It  is  assumed  that  the  vehicle  is  mounted  to  rotate  about  its  center 
of  gravity,  which  is  centered  in  the  design  centerline  of  the  solar  bank, 
and  the  vehicle  is  axially  symmetric.  Thus  the  outline  of  the  vehicle 
projected  on  the  solar  bank  is  single  valued  in  any  quadrant  and  is  sym¬ 
metrical  about  the  vertical  centerline  on  the  solar  bank. 

The  program  projects  the  vehicle  outline  and  tests  according  to  the 
option  to  determine  if  the  solar  module  should  be  on. 

The  geometry  of  the  vehicle  is  constant  for  any  run.  The  pitch 
angle  is  variable  and  is  orbit  and  orientation  dependent.  The  eclipse 
condition,  all  solar  modules  off,  is  orbit  dependent  only.  Vehicle  roll 
does  not  enter  these  calculations  because  axial  symmetry  is  assumed. 


4.0  SECONDARY  ILLUMINATION 


We  have  considered  in  previous  sections  the  means  of  calculating 
orbital  motion  and  the  main  simulation  variables,  roll  and  pitch,  which 
properly  orient  the  vehicle  to  the  solar  bank.  We  must  now  look  into 
the  secondary  illumination  fields,  albedo  and  earth  radiance. 

Roll  of  the  vehicle  takes  place  within  the  simulation  cell  to  keep  the 
vehicle  "down"  direction  aimed  at  earth's  center*.  Since  roll  is  relative 


Move  strictly,  to  keep  a  particular  reference  plane,  containing  the 
roll  axis,  aligned  with  the  center  of  rotation  of  the  body. 
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to  the  "birdcage"  array,  this  means  that  the  down  direction  moves  with¬ 
in  the  birdcage,  and  this  motion  is  just  the  roll. 

As  we  will  bring  out  in  what  follows,  the  real  secondary  field  has 
an  extended  area  source,  the  earth's  surface  or  part  of  it;  hence,  the 
flux  arriving  at  the  vehicle  in  orbit  is  a  continuous  function  of  angle  at 
the  vehicle.  But  the  birdcage  is  a  system  of  six  discrete  broken  line 
arrays,  equally  spaced  about  the  vehicle.  Thus  the  best  that  can  be 
hoped  for  in  simulation  is  to  match  integrals  of  flux.  That  is,  by  sub¬ 
dividing  the  vehicle  into  six  zones,  it  should  be  possible  to  set  the  six 
birdcage  lines  at  power  levels  such  that  the  integrated  flux  into  a  zone 
on  the  vehicle  matches  the  integrated  flux  into  the  same  zone  on  an 
orbiting  vehicle  exposed  to  the  real  earth  illumination  source.  This  is 
the  principle  problem  of  albedo-radiance  simulation.  Before  such 
matching  can  be  accomplished,  it  is  necessary  to  know  the  natures  of 
the  real  albedo  and  radiance  functions  and  of  the  simulation  element 
distribution  function. 

In  our  approach  to  this  problem,  we  assume  the  vehicle  is  axially 
symmetric  and  tangent  oriented.  Then  we  will  develop  profiles  of  flux 
incident  on  a  band  around  the  girth  of  the  vehicle.  The  profile  will 
depend  on  position  in  orbit  and  the  position  on  the  vehicle.  So  we  ob¬ 
tain  typically  a  single  profile  for  each  pair  of  9 a  ,  values.  We  then 

assume  that  the  six  birdcage  element  patterns  are  of  identical  form, 
differing  only  by  equal  phase  shifts  in  their  arguments  and  by  their  am¬ 
plitudes.  The  six  superimposed  birdcage  patterns  also  are  shifted  rela¬ 
tive  to  the  profile  by  an  amount  equal  to  the  roll. 

The  birdcage  element  illumination  pattern  at  the  vehicle  depends  on 
the  relative  size  of  the  birdcage  and  vehicle  and  on  any  optical  surfaces 
involved.  This  is  not  known  in  sufficient  detail  at  this  time  but  will  be 
assumed  to  have  the  desirable  properties  of  uniformity  over  one-sixth 
the  girth  of  the  vehicle  with  small  "tails"  extending  into  neighboring 
zones. 

There  is  also  an  axial  distribution  pattern  inherent  to  the  design  of 
heater  elements  in  essentially  linear  arrangement.  If  all  elements  were 
equally  excited  along  the  line,  a  maximum  flux  occurs  opposite  the  mid¬ 
point  of  the  array.  Thus  to  even  out  the  flux  at  a  fixed  distance  from 
the  source  requires  the  elements  be  excited  roughly  in  proportion  to  the 
reciprocal  of  the  pattern  intensity  for  a  uniformly  excited  array.  The 
end  elements  will  thus  be  hotter  than  the  middle  elements.  This  will  be 
treated  in  detail  in  section  4.  1. 
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4.1  SECONDARY  ILLUMINAT ION  SIMULATION 

In  this  section,  we  will  examine  the  radiation  field  produced  by  an 
array  of  linear  elements  to  obtain  some  idea  of  the  longitudinal  distribu¬ 
tion  of  radiation  produced  by  the  "birdcage. "  We  confine  attention  to 
the  cylindrical  portion  of  the  array  about  a  cylindrical  model.  We  will 
assume  that  the  longitudinal  and  lateral  distributions  are  independent. 

We  look  first  at  the  longitudinal  distribution. 


Now  consider  an  infinitely  long  line  source  with  flux  w  emitted  per 
unit  length  41  (Fig.  16).  An  element  of  area  dA  at  distance  R  is  oriented 
so  that  it  is  parallel  to  the  line  source  and  its  normal  intercepts  the  line. 
The  flux  dW  incident  on  dA  is 


dW 


4-  oc 


Wd  A 


f 


d  di 


The  polar  coordinate  description  of  a  line  is 
hence, 


d  W 
d  A 


it/, 2 

f  cos 2  Odd 
J~n/2 


ttVI 
2  R 


p  cos  6  =  R 


=  I 


oo 


and  t  ?=  R  tan,0, 


the  intensity  of  flux  for  an  infinitely  long  line. 
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For  a  finite  line,  the  previous  integral  is  evaluated  between  finite 
limits : 


w_ 

R 


6 


f. 


cos2  did  =  If 


And  for  an  array  of  N  linear  line  segments  (i  =  1,  2,  • 


w 

R 


If  ,  n 


N)  ,  we  have 


V  A 

The  limits  of  integration  $i,  0 i  are  determined  by  the  location  of 
the  area  element  dA  relative  to  the  ends  of  the  i th  element.  It  is  con¬ 
venient  to  normalize  relative  to  the  infinite  line  flux  1  .  Then  the  inte¬ 
grals  are  of  form 


If  N  elements,  each  of  length  L,  spaced  S  on  centers,  are  con¬ 
sidered,  we  obtain 


jf  ,  N 
i 

03 


R  L 


R  +  ( 


—  X  )  (  X; 


-x) 


R  (  X}  —  x  ) 

2  .  A  2 

R  +  (  XJ  —  X  ) 


R  ( x*  —  x  ) 

2  V  2 

R  +  (  X;  -  X  )  _ 


where 


A  V  A  A  ,, 

xj  —  Xj  +  JLj  ,  x i  4-1  —  xj  +  o 


The  results  for  various  values  of  R  are  plotted  in  Fig.  17  for  ten 
elements  having  L  =  0.9,  S  =  1  .  For  comparison,  the  results  for 
R  =  10,  L  =  10,  S  =  10  are  plotted,  showing  that  at  a  little  distance,  the 
small  discontinuities  in  an  interrupted  line  source  are  washed  out,  and 
the  flux  is  essentially  that  of  a  single  source  of  slightly  less  strength  of 
the  same  overall  length.  (See  Ref.  6  for  results  of  additional  work  at 
AEDC  on  the  line  source  arrays. ) 

The  important  feature  here  is  the  gross  shape  of  the  curves,  maxi¬ 
mum  at  center  and  dropping  at  ends.  This  means  that  the  inverse  of 
this  function,  or  of  its  piecewise  integrals,  will  be  needed  to  weight  the 
intensities  of  the  individual  elements  of  the  source  to  provide  a  uniform 
illumination  region  at  R  distance  from  the  source. 
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Now  the  actual  birdcage  elements  are  not  simple  linear  sources; 
they  are  rods  of  definite  size  at  the  foci  of  cylindrical  reflections  of 
roughly  piecewise  parabolic  segment  cross  section.  The  reflector 
design  is  aimed  at  defining  rather  sharply  the  lateral  illumination  on  a 
specific  body  (nominally  a  right  circular  cylinder  plus  cone).  The 
reflectors  thus  have  the  principle  characteristic  of  cylinder  optical 
elements;  they  are  astigmatic.  This  means  that  any  point  in  object  space 
(on  the  source  element)  will  be  imaged  (real  or  virtual)  as  a  line  in 
image  space.  This  fact  adds  to  the  validity  of  the  foregoing  discussion. 

The  real  line  source  elements  are  separated  by  plane  baffles  per¬ 
pendicular  to  the  line,  perhaps  reflectors,  to  limit  the  angular  extent 
of  radiation.  Since  their  dimensions  are  small  relative  to  the  length  of 
line  segment,  they  effect  the  distribution  only  slightly,  probably  insig¬ 
nificantly.  They  most  strongly  influence  the  flux  arriving  at  a  given 
point  from  most  remote  elements,  which  contribute  very  little  in  any 
case. 

The  result  of  baffling  is  indicated  in  Fig.  18  for  a  single  element. 
The  main  features  are  slight  increase  and  flattening  of  peak,  an  addition 
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to  shoulder,  a  reduction  of  heel,  and  a  finite  termination  of  the  tail  of 
the  distribution,  all  minor.  The  N  element  distribution  can  be  obtained 
by  superposition. 


Source  Element 


Limited  to  Single  Reflections  by  Perfect  Mirrors;  Subsequent  Reflections 
Contribute  Much  Less. 

Fig.  18  Approximate  Effect  of  Reflecting  End  Baffles  on  Longitudinal  Flux  Distribution  for  a 
Single  Line  Source  Element 


The  source  elements  are  in  reality  small  rods  rather  than  lines. 

We  have  used  the  notation  W  for  intensity  per  unit  length  from  a  line. 

For  a  directionally  homogeneous  radiating  line,  whatever  radiation 
leaves  the  line  also  passes  through  a  cylinder  of  sufficiently  small  diam¬ 
eter  concentric  with  the  line.  If  the  rod  dimensions  are  small  relative 
to  other  lengths,  the  assumption  of  ideal  line  source  is  a  good  approxi¬ 
mation. 

We  have  assumed  that  the  lateral  and  longitudinal  distributions  are 
independent.  This  is  an  essential  point  with  respect  to  computability. 

If  it  is  necessary  to  control  all  elements  of  the  birdcage  independently, 
we  shall  have  to  solve  frequently  a  set  of  192  equations  in  192  unknowns 
to  find  the  amplitudes  -at  which  to  excite  the  192  birdcage  elements.  If 
the  assumed  independence  is  valid,  this  reduces  to  6  equations  in  6  un¬ 
knowns  for  the  lateral  distribution  and  to  the  application  of  a  previously 
evaluated  weighting  function  for  the  longitudinal  distribution. 
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We  must  mention  some  design  problems  posed  by  the  birdcage  array. 
First,  a  vehicle's  base  will  intercept  secondary  illumination,  but  base 
illumination  does  not  occur  in  the  birdcage  design,  except  incidentally. 

Second,  a  real  vehicle  of  say  cylinder-cone  configuration,  when 
about  horizontal,  will  have  a  zone  "on  top"  the  cylinder  which  cannot 
see  the  earth;  whether  the  top  of  the  cone  has  such  a  zone  depends  on 
the  altitude  and  cone  angle.  If  the  half  cone  angle  is  greater  than  the 
horizon  depression  with  altitude,  the  cone  will  be  illuminated  when  the 
cylinder  is  not.  But  the  birdcage  elements  will  not  discriminate;  if  an 
element  is  "on"  to  illuminate  the  cone  it  will  also  illuminate  the  neigh¬ 
boring  cylinder  which  should  not  be  illuminated.  If  the  vehicle  nose  is 
ogival  or  another  smoothly  faired  shape  having  no  discontinuity  of  slope 
where  it  joins  the  cylinder  body,  this  is  not  so  severe  a  problem. 

Third,  the  sharpest  defined  boundary  for  illumination  with  line 
sources  and  cylindrical  reflectors  are  the  lateral  boundaries,  parallel 
to  line  arrays.  All  other  boundaries  are  smeared  by  the  astigmatic 
reflectors.  Thus  it  is  not  possible  to  sharply  define  the  illumination  on 
a  pointed  nose  piece.  A  practical  solution  here  would  be  to  define  the 
nose  by  means  of  a  cone  frustrum  roughly  matching  the  nose  shape.  The 
nose  illumination  pattern  will  be  inefficient  in  any  case,  and  consider¬ 
able  illumination  will  pass  directly  to  the  test  chamber  walls  and  increase 
the  heat  load. 


4.2  “BIRDCAGE"  PATTERNS  AND  MATCHING 

Assuming  that  we  know  the  combined  illumination  profile  on  a  vehicle 
in  orbit,  we  must  go  through  some  computations  to  find  the  birdcage 
excitation  required  to  effectively  reproduce  this  profile  on  the  test  vehicle 
(Figs.  19,  20,  21,  and  22). 

If  there  are  M  discrete  sources  of  amplitudes  a;  (i  =  1,  2,  ....  m)  , 

having  distribution  functions  (from  source  i  into  region  k  on  vehicle) 
and  the  profile  function  fk  (from  real  source  into  region  k),  the  best  we 
can  hope  to  do  is  to  match  the  total  input  into  each  of  the  M  regions,  by 
proper  selection  of  t  matching  real  input  to  simulating  input  over 
each  region  k  .  We  are  limited  to  M  regions  on  the  vehicle  by  the  M  dis¬ 
crete  sources,  assumed  so  arranged  as  to  cover  the  vehicle. 

We  require  for  the  area  Sk  of  each  k  region,  k  =  1,  2,  •  •  •  M, 

/  fkdS  =  L  ai  /  Sik  dS 
sk  1  =  1  sk 
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Fig.  20  Profile  Matching  with  Integrals 
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Relative  Orbit,  Earth  about  Vehicle 


Sun 


"Down"  Points 
to  Earth  Center 


Birdcage 
Fixed  Relative 
un 

Fixed 
Hexagon 
Approx, 
to  Body 


Since  "Down"  is  Earth  Fixed  in  Vehicle 
and  Hexagon  Pattern  is  Sun  Fixed,  Relative 
Motion  (Roll)  Occurs  between  Hexagon 
Pattern  and  Vehicle. 


Fig.  21  Six-Point  Profile  Matching 
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Fig.  22  Profile  Matching  at  Six  Points  Only 
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The  git  are  fixed  by  birdcage  element-reflector  geometry.  The  fy 
are  orbit-orientation  functions  and  depend  on  position  on  the  vehicle. 
They  depend  on  location  on  the  vehicle  surface  and  on  the  time,  through 
orbit  position.  We  require  a  set  of  a;  for  each  simulation  interval 
(nominally  two  minutes  real  time). 

The  direct  solution  of  this  system  of  M  equations  in  M  unknowns 
leads  to  appearance  of  negative  aj’s,  obviously  meaningless.  By  a 
"least  squares"  approach,  we  obtain  the  system  of  normal  equations 

M  M  M 

L  ai  I  Gik  Gjk  =  Z  Fk  Gjk  >  j  =  1 ,  2 ,  .  .  .  M 

i=  1  k  =  1  k  =  1 

where  we  have  used  the  definitions 

Gq  =  f  g..  dS  Fk  =  /  fk  dS 

^k  Sk 

In  vector -matrix  notation,  we  have 

{  S  }  a  =  V  or  a  =  !  S  }  V 
where  IS}  is  a  matrix  of  elements 

S i j  =  2  G ij  Gik 

—  1 

[  S  }  is  its  inverse 

V  is  a  vector  of  components  Vj  =  2  Fk  Gjk 

a  is  a  vector  of  components  aj 

Ordinarily  the  matrix  elements  are  functions  of  geometry  and  orien¬ 
tation  of  the  vehicle -discrete -source  system,  and  the  Vj  are  functions 
of  vehicle  position,  orientation,  altitude,  etc.  ,  as  well.  We  may  say 
the  Vj  are  orbit  dependent  (through  the  fk  function)  and  roll  dependent 
(through  the  g;k  function)  since  they  depend  on  roll  of  vehicle  within  the 
birdcage. 

But  if  the  vehicle  is  perfectly  axially  symmetric  in  shape  and  sur¬ 
face  properties,  the  g;k  become  independent  of  roll;  hence,  the  matrices 
}  S  }  and  {  S  }  have  constant  elements,  and  the  Vj  are  independent  of  roll. 
This  point  is  of  considerable  significance. 

If  the  S  matrix  is  constant,  we  may  immediately  invert  it  one  time 
for  j  S }  and  program  the  evaluation  of  M  equations  in  M  unknowns  at 
each  point  in  orbit.  But  if  the  iSl  matrix  is  variable,  we  must  invert 
it  for  { S }  at  each  point  in  orbit,  then  solve  the  M  equations  in  M  un¬ 
knowns. 
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For  the  most  rapid  computing,  if  the  resulting  errors  are  not  ex¬ 
cessive  in  practice,  we  may  replace  the  circular  section,  on  which  the 
profiles  are  based,  by  a  hexagonal  section.  Each  leg  of  the  hexagon 
directly  faces  one  of  the  birdcage  line  elements,  and  we  ignore  any 
"tails"  of  the  element  distribution  function.  Then  the  amplitudes  ai  are 
directly  associated  with  six  corresponding  points  on  the  profile. 

Since  the  profile  pattern  changes  continually  with  orbit  position, 
and  roll  changes  similarly,  we  must  still  evaluate  six  new  profile  am¬ 
plitudes  at  each  point  in  orbit.  These  are  unique  calculations  at  each 
point  in  realistic  orbits  and  repeat  periodically  only  for  a  static  universe 
orbit  model  in  which  the  orbit  period  is  an  exact  multiple  of  the  simula¬ 
tion  interval  (2  min). 


4.3  ILLUMINATION  INTEGRALS 

To  compute  earth  Albedo  and  radiance  integrals  for  a  surface  ele¬ 
ment  having  arbitrary  orientation  and  position,  we  make  the  following 
assumptions: 

1.  that  Albedo  is  a  uniform  property  of  earth's  surface, 

2.  that  the  earth  is  a  sphere, 

3.  that  the  earth's  surface  is  diffusely  reflecting,  and 

4.  that  the  earth  has  no  atmosphere. 

Is  Solar  constant,  intensity  of  solar  radiation  at  earth 

ae  Albedo,  fraction  of  solar  constant  reflected,  a 
surface  property 

The  solar  radiation  incident  on  an  area  element  AAt  having  its  normal 
inclined  at  angle  6e  to  sun  direction  Is 

f  Is  cos  0e  A  A.  for  cos  0e  >  0 

AIe  =  I 

l  0  for  cos  0e  <  0 

Of  this  a  fraction  ae  is  reflected  diffusely  by  AAt ;  hence,  the  in¬ 
tensity  per  unit  solid  angle  AI^,  in  a  direction  inclined  at  angle  to  the 
surface  normal  is 

AI^  =  a^—  cos  0e  COS  iff  AAj. 

The  intensity  included  in  solid  angle  Aoj  is 

AI^j  =  a^s —  cos  9e  cos  i ft  AAt  A  a 
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An  area  element  AA2  J  at  distance  Pe  from  AAt  ,  having  its  normal 
inclined  at  angle  <f-  to  the  direction  of  Pe  >  intercepts  a  solid  angle 
(Fig.  23). 

.  A  A,  cos  £ 

Aco  =  - - — ~ - - — 

Pe 


Hence  the  intensity  arriving  at  AA2  is 


AI2 


is 


COS  6e  COS  xjj  COS  £  AAi  A  A2 


Fig.  23  Solid  Angle  Geometry  for  Albedo  and  Earth  Radiance  Calculation 


Let  the  area  element  AA2  be  located  at  altitude  h  above  earth  of 
radius  Re  (Fig.  24).  From  this  point,  the  portion  of  the  earth  that  can 
be  seen  is  confined  within  a  horizon  circle.  The  angle  y  between  the 
direction  of  earth  center  and  the  horizon  circle  is  defined  by 

sin  y=Re/(Re+h)  (0  <  y  <  ) 

At  the  earth  center,  let  6S  be  the  angle  between  the  direction  to 
Aa2  and  the  sun  direction;  let  6e  be  the  angle  between  the  area  element 
AA!  on  earth  and  the  sun;  let  4>e  be  the  angle  between  the  directions  of 
AAj  and  AA2 . 
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siny  =Re/(Re  +  h) 

cos  0e  =  cos  0  cos  0S  +  sin  0  sin  0S  cos  a 
sin  ip  =  sin  p/sin  y 

i|i  =p  +  0 

Fig.  24  Coordinate  Relations  for  Albedo  and  Earth  Radiance  Calculation 

At  the  element  AA2 ,  xjj  is  the  angle  between  the  normal  to  At  and  the 
direction  of  AA2  as  before. 

At  the  element  AA2  ,  let  /3  be  the  angle  between  the  direction  to  earth 
center  and  AA2 .  Angle  /3  is  a  "nadir"  angle. 

Then  the  following  relations  hold: 

'A  =  /3  +  <f> 

sin  ip  =  sin  j8/ sin  y 

cos  Qe  =  cos  0§  cos  <f>e  +  sin  @s  sin  <Ae  cos  a 

where  a  is  the  angle  about  the  line  from  AA2  to  earth  center  measured 
from  the  plane  including  this  line  and  the  sun.  Angle  a  is  the  azimuth 
angle  (Fig.  25). 

Now  the  element  is  described  in  spherical  coordinates  having 
polar  axis  along  the  earth-to-AA2  line,  longitude  a ,  and  co -latitude  <Ae: 

AAX  =  Re2  sin  <Ae  d  a  d  pe 
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l  he  a,  p  Coordinate  System  at  dA2 


P  =0 

cos  i  =cosPcosPN  +  sinPsinPNcos(a-aN) 


Fig.  25  Albedo-Radiance  Integration  Coordinates 


Similarly,  we  may  describe  a  spherical  area  element  in  terms  of 
a  ,p  ,  and  f>e  from  AA2  : 

AAp  =  pe2  sin  /3  d  j3  d  a 


Hence  any  point  may  be  described  by  either  (Re,  a)  or  (pe  ,  (3  ,  a) 
at  AAj  ,  and  we  find  that  AAp  and  AAt  are  related  by  simple  projective 
properties: 

A At  cos  i ft  =  AAp 

so  that  A  A,  =  ^  sin^  da 

COS  XjJ 


We  may  now  write  the  intensity  at  AA2  caused  by  reflection  from  AAt 


Ai2 
A  a2 


COS  6e  COS  I jj 


^  Sin-^..  da  d  (3 


ae  Is 
77 


cos  de  cjjs  g  sin  /3  d  j3  d  a 
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Integrating  over  the  interior  of  the  horizon  circle,  we  obtain 

T  t  ax 

.  2  -  (A  A2 ,  0S  >  y  >  0  =  3e-  S —  J  I  cos  (d,  /3,  &S,  y)  cos  £  sin  /3  d/3  da 

AA*  77  £  =  o 

where  amin  ,  amax  may  be  functions  of  /3 ,  and  there  may  be  several  dis¬ 
tinct  regions  or  functions. 

In  the  spherical  coordinates  a  ,  j3  ,  any  orientation  of  surface  ele¬ 
ment  may  be  described  by  the  direction  angles  of  its  normal,  an  ,  /3n 
(Fig.  26).  Then  the  angle  £  is  found  from 

cos  f  =  cos  /3  cos  /3n  +  sin  /3  sin  /3n  cos  (a  —  apf ) 

Since  0  =  ,/r  -  / 3,  and  sin  ifi  =  sin  /3/ sin  y 

cos  0e  =  sin  i/f  (cos  0s  sin  /3  +  sin  0$  cos  /3  cos  a) 

+  cos  i /r  (cos  0s  cos  /3  —  sin  0s  sin  /3  cos  a) 

The  intensity  integrand  is  completely  expressible  in  the  two  vari¬ 
ables  a  ,  j3  ,  and  the  configuration  parameters  0s  ,  y  ,  a^  ,  /3N  . 

At  this  stage  it  is  possible  to  integrate  numerically  by  letting  a 
range  from  0  to  2n  and  /3  range  from  0  to  y ,  provided  that 

cos  £  >  0 

cos  0e  >  0 

COS  if'  >  0 

and  using  (=0)  for  any  (a,  /3)  violating  these  conditions.  This  will  be  dis¬ 
cussed  later  in  Appendix  IV. 


Normal,  Cone 


Fig.  26  Area  Elements  for  Cylinder-Cone 
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Ignoring  for  the  moment  the  constant  ae  Is/ 77,  we  have  the  following 
terms  to  be  integrated  over  a  and  /3  : 


1. 

COS  /S ]\J 

COS 

Os 

cos2  /3 

sin  /3 

COS  1 jj 

2. 

—  cos  /3n 

sin 

Os 

cos  /3 

sin2  /3 

cos  1/7 

COS 

a 

3. 

cos  /3]\j 

cos 

Os 

cos  /3 

sin2  /3 

sin  ip 

4. 

cos  /3j\[ 

sin 

Os 

cos2  /3 

sin  /3 

sin  ip 

cos 

a 

5. 

sin  /Sn 

cos 

Os 

cos  aN 

sin2  /3 

cos  j3 

cos  ip 

cos 

a 

6. 

sin  /S jvf 

cos 

Os 

sin  aft 

sin2  /3 

cos  /3 

cos  ip 

sin 

a 

7.  -  sin  /3n  sin  cos  «N  sin3  /3  cos  ip  cos2  a 


8. 

—  sin 

0N 

sin  0s 

sin  aft 

sin3  / 8 

cos 

<A 

sin 

a  cos  a 

9. 

sin 

cos  0s 

cos  aN 

sin3  /S 

sin 

-A 

cos 

a 

10. 

sin 

/3n 

cos  0s 

sin  gen 

sin3  /3 

sin 

-A 

sin 

a 

11. 

sin 

sin  0 s 

cos  a]\j 

sin2  /3  cos 

0 

sin 

'A 

cos2 

a 

12. 

sin 

/3  N 

sin  0s 

sin  a]\f 

sin2  /3  cos 

£ 

sin 

-A 

sin 

a  cos  a 

Note: 

sin  ^ 

=  sin  /3  /  sin 

y 

As  long  as  there  are  no  boundaries  of  earth  surface  for  which 
a  =  a  { j3  )  ,  so  that  a  can  range  from  0  to  2n  ,  we  may  integrate  relative 
to  a  and  obtain  simple  results.  Integrals  (1,  3)  do  not  contain  a  so  the 
integration  results  in  a  factor  2 n .  Integrals  (2,  4,  5,  6,  9,  10)  contain 
only  sin  a  or  cos  a  and  vanish.  Integrals  (7,  8,  11,  12)  contain  cos2  a  or 
cos  a  sin  a  and  result  in  a  factor  of  n  . 

These  conditions  are  satisfied  as  long  as  we  have  both 

*  -f-  -  y 

Os  ^  y 

If  /3n  ^  n/2  +  y  or  %  >  77  -  y  ,  the  earlier  conditions  on  cos  £  or 
cos  0e  are  violated  and  the  entire  integral  (I2/ AA2)  vanishes. 

For  -f - y  <  /3n  <  -f -  +  Y 

and/or  y  <  0s  <  n  -  y 

there  exist  boundaries  of  form  a(/3),  and  the  integration  becomes  com¬ 
plicated.  The  integration  is  bounded  by  arcs  of  one,  two,  or  three 
curves  of  a{/3),  whose  intersections  are  generally  given  by  implicit 
functions.  A  first  integration  may  be  done  formally;  expressions  result 
for  which  the  integrals  are  not  available  in  closed  form. 
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Numerical  integration  may  be  accomplished  in  an  easily  compre¬ 
hended  manner  by  referring  to  the  earlier  integral  expression.  The 
product  cos  f  cos  9e  sin  j3  may  be  calculated  term  by  term,  and  in  addi¬ 
tion,  the  expression  cos  i fj  can  be  evaluated  to  ensure  that  the  conditions 


are  satisfied. 


COS  f 

cos  9e 

COS  I jj 


0 


For  some  value  combinations  of  0$  the  integrations  can  be 

carried  out.  Five  special  cases  are  given  in  Appendix  II.  A  typical 
body  for  which  profiles  are  easily  obtained,  assuming  circular  orbit, 
tangent  orientation,  is  shown  in  Fig.  26. 


4.4  RADIANCE  FUNCTIONS 
4.4.1  Earth  Radiance  Integrals 


The  earth  is  taken  to  be  a  warm  body,  uniform,  having  a  flux  Ir 
issuing  from  each  element  of  area.  The  corresponding  flux  arriving  at 
the  element  of  area  AA2  is  obtainable  from  the  albedo  integral  by  setting 
cos  0e  =  1  and  replacing  constant  Is  by  Ir  . 


la 

A  A2 


cos  f  sin  /3  d  /3  da 


and  cos  f  =  cos  /3  cos  /3n  +  sin  /3  sin  /3]\j  cos  (a  -  a^)  ,  as  before. 


The  integration  is  given  in  Appendix  III  with  changes  in  form  depend¬ 
ing  on  the  boundary  relationships.  The  closed  forms  are  available  because 
of  the  great  simplification  of  integration  over  a  region  bounded  by  one  or 
two  arcs  and  the  simpler  integrand. 

4.4.2  Combined  Secondary  Illumination 


The  illumination  intensity  to  be  simulated  by  the  birdcage  is  the  com¬ 
bined  Albedo-earth  radiance  functions 


]§_ 

A  A2 


IS  ae 
77 


COS  f  sin  j3 


d/3 


d  a 


Ir  J"  J  cog  ^  gjn  p  J  jg  d  a 
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A  factor  (Eff)  is  included  in  the  second  term  for  purposes  of  correct¬ 
ing  the  spectral  differences  between  earth  and  the  incandescent  heater 
elements  of  the  birdcage.  This  factor  may  be  in  fact  a  function  of  posi¬ 
tion  on  the  vehicle,  hence  a  function  of  the  roll  angle.  It  is  a  surface 
property  of  the  vehicle. 


5.0  CONCLUSIONS 


Provided  that  some  simplifications  are  allowed,  it  is  possible  to 
organize  a  digital  computer  system  of  7074  programs  which  will  provide 
realistic  orbit  calculations  and  control  of  all  the  Mark  I  simulation 
parameters  in  a  small  fraction  of  real  time.  Any  violation  of  the  sim¬ 
plifying  assumptions  will  require  some  special  programming  and 
generally  more  computer  time.  Most  of  the  program  types  have  been 
written  and  run  in  connection  with  Mark  I  projects  but  only  piecemeal. 

A  logical,  next  step  in  this  sequence  is  the  organization  of  a  compre¬ 
hensive  program  to  do  all  required  calculation  for  really  typical  test 
vehicle-orbit  configurations  with  as  high  realism  as  computer-time-to- 
real-time  ratios  will  allow.  Perhaps  a  more  immediately  practical 
approach  would  be  to  develop  a  fast  computing  program  system  using 
rather  gross  approximations  and  to  develop  a  slower,  more  accurate, 
highly  realistic  study  system  for  evaluating  the  overall  results  of 
approximations  on  the  quality  of  Mark  I  performance  and  test  results. 
Having  then  the  two  extremes  programs  available,  future  efforts  would 
be  aimed  at  obtaining  an  efficient  practical  compromise. 
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APPENDIX  I 

A  LIST  OF  IBM  7074  PROGRAMS  CONCERN  ED  WITH  MARK  I  ILLUMINATION  SIMULATION 


Program  No. 

Title 

05027 

Interplanetary  Trajectory  Study  Mark  I 

05028 

Orbit  Calculations 

05032 

Interplanetary,  Trajectory  Study  Mark  I  (RSW-1) 

05036 

Double  Precision  Orbits 

05036 

Orbit  Subroutine 

05045 

Encke  Orbit 

05048 

Interplanetary  Trajectory  (RSW-2) 

05049 

Interplanetary  Trajectory  (RSW-3) 

05050 

Orbit  III 

05051 

Radiation  Study 

05064 

Albedo  Ray  Trace  (1) 

05073 

Albedo  Profiles 

05087 

Total  Heat  Input,  Cone  or  Cylinder 

05090 

Themis  Translation 

05092 

Roll  Pitch  Plot 

05109 

Albedo  Plot 

05122 

Albedo  Profile  Transformation 

05123 

Albedo  Integration 

05133 

Albedo  Integral,  Transition 

05205 

Solar  Simulation  Module  Honeycomb 

05064 

Albedo  Ray  Trace 

05212 

Albedo  Ray  Trace  Revision 

05216 

Albedo  JMirror  Design 
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APPENDIX  II 

SPECIAL  CASES  OF  ALBEDO  INTEGRAL 


CASE  1 

AA2  horizontal  (/3n  =  0)  and 
located  on  line  between  earth  and 
sun  (Os  =  0)  . 

This  gives  the  maximum  value 
of  the  integral  and  is  very  useful 
for  order  of  magnitude  estimates. 

With  /3n  =  0  >  Os  =  0  ,  integrands 
(1,  3)  are  the  only  contributors: 

1(a)  2  tt  cos  /3n  cos  0g  cos2  /3  sin  /3  cos  ip 
3(a)  cos  j8n  cos  0g  cos  /3  sin2  /3  sin  i jj 

become 

1(b)  2  tt  cos2  /3  sin  /3  1  —  S*n2  ^ —  d/3 

3(b)  2  tt  cos  jS  — ^-r — —  d/3 

sin  y 

which  are  simply  transformed  to 

1(c)  2  tt  — s-^  — V  cos2  /3  —  cos2  y  d  (  —  cos  /3  ) 

3  (c)  2  77  ~s?n  y  d  Sin  ^ 

3(c)  integrates  simply  to  -y  sin3  y  over  the  grange  (0,  y)  while  1(c) 
integrates  to 


dyS 

d/S 


77  i 

r  .3 

1  sin  y  + 

cos2  y  sin  y 

cos4  y  / 

'  cosy  \ 

~T~  ] 

2  + 

2  ln  ( 

v  1  +  sin  y y 

So  for  the  "hot  spot"  (/3  =  0,  %  =  0) 


A  A, 


=  ae  IS 


.  3  COS  y  .  2  7 

sin  y  +  - — ! —  <  sin  y  +  cos  y  in 


( 


cos  y 
1  +  sin  y 
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CASE  II 

AA2  is  located  over  the  termina¬ 
tor  (0s  =  n'/2 )  and  horizontal.  Inte¬ 
grands  (2,  4)  survive.  The  range  of  a 
is  (  —  77/2 ,  +  n/2 )  . 

2.  -  cos  /3  sin2  /3  cos  ifi  cos  a  d  a  d /3 

4.  cos2  sin  (3  sin  if/  cos  a  da  d f3 


cos  /3  sin2  /3  / 1 


2  (a) 

=  -  sin2  /3  xj 1  - 

integrating  to1 


|3 


sin  y 


d  ( sin  a)  d  /3 


/3 


sin  y 


d  (sin  a)  d  (sin  /3 ) 


Sun 


2  (b)  - 


sin  y 


-  sin/3  (sin2y  -  sin2 /3)/2  +  —  \  sin  (3  (sin2  y  -  sin2  j8)/a 


+  sin  y  sin 


-1  (- 
\  si 


/3 


sin  y 


sin  y 


4(a) 


0  „  sin2  /3  ,  ,  .  %  .  ^ 
3*  P  sin  v  -  d  ^sln  d  P 


sin  y 
1 

2  sin  y 
1 


4  sin  y 

for  this  case,  /6n  =  °>  %  =  n/2 


J8  -  4-  1/3  + 


y  —  sin  y  cos  y  (cos2  y  —  sin2  y) 


A  A, 


ae  *S 
4  77 


1 


sin  y 


y  -  sin  y  cos  y  +  cos  y  sin  y 


sin  y 
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CASE  III 


6s  =  0,  jSn  =  77-/2  .  Integrands 
5,  6 ,  9,  10  remain. 

5.  cos  otN  sin2  P  cos  J 0  cos  'A  cos  a 

6.  sin  ajsf  sin2  /3  cos  /3  cos  ifj  sin  a 

7.  cos  a]\[  sin3  /3  sin  ip  cos  a 

8.  sin  aj\f  sin3  /3  sin  i p  sin  a 


which  combine  as 

(5,  6)  (a)  sin2  6  -7- ^  V  sin2  y  -  sin2  B  cos  (a  -  ajvj ) 

sin  y 


(7,  8)  (a) 


P 


sin  y 


cos  (a  —  ) 


and  (a  -  aN)  ranges  (-  u/2,  +  n/2)  integrating  to  a  factor  2. 


Then  with  respect  to  ft  the  integration  gives,  similar  to  case  II, 
(5,  6)  (b)  +  sin3  y 

and  — -  -7-  y - i-  sin  y  cos  y  +  (sin  y  cos2  y  —  cos  y  sin2  y  ) 

sin  y  L  8  2  '  '  J 

SO 


Is 


A  A, 


sin  y  + 


1 


7T  sin  y 


(  3y  —  4  sin  y  cos  y  +  sin  y  cos2  y  —  cos  y  sin2  y) 


CASE  IV 


#S  =  it/2,  £n  =  tt/2  ,  ctN  =  0  . 
Integrands  7  and  11  remain  with 
a  ranging  -n/2,  +  n/2 

(7)  -sin 3 /3  cos  ^  cos2  a  da  d /3 

(11)  sin2  jS  cos  /3  sin  i/'  cos2  a  da  d /3 

becoming,  through  a. 


sin2  /3 

s  in2  y 


sin3  /3 
sin  y 
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and,  through  /3  , 


CASE  V 

%  =  u/2,  j3 n  =  n/2,  on  =  n/2  .  a  ranges 
(0,  77/2)  giving  77/ 4  in  a  .  Integrands  8  and  12 
survive. 

(8)  -  sin3  B  Jl  —  - 5-^ —  sin  a  da  dl3 

v  sin  y 

( 12  )  sin2  (3  cos  (3  sin  xft  sin  a  cos  a  d  a  d  j3 

(12)  goes  simply  to 

(12)  (a)  +  Sin40  7 - 

8 

(8)  gives 

(8)  (a)  — — 7- -  -  sin  y  (  +  _.cos  \  +  cos2  y  (  1 

so  [  12  (a)  +  8  (a)  ]  gives  3e — 

S  AA2  /  77 
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APPENDIX  III 

EARTH  RADIANCE  INTEGRAL 


We  must  evaluate 

XP  max  /»Aamax 

J  cos  £  (a,  (3  ,  aN  ,  /3n  )  sin  jS  d/3  d  A  a 

fJmin  ^amin 

For  all  [  a,  jS  ]  for  which,  with  A  a  =  a  -  , 

cos  £  =  cos  /3  cos  jSn  +  sin  /3  sin  /3]sf  cos  A  a  >  0 

Integrating  with  respect  to  A  a 


(1) 


(2) 


fimax 


I 


Ja  m  m  a  x 

[ 


cos  /3]\j  cos  /3  sin  /3  A  a  +  sin  j8j([  sin2  /3  sin  A  a 


A  a„ 
A  am 


d/3  (3) 


By  symmetry  and  oddness  of  Aa,  sin  A  a  ( both  zero  at  Aa  =  0) 


P, 


I  = 


*  L  [ 


cos  /3n  cos  ^3  sin/3Aamax  +  sin sin2 /3  sin  Aamax  d /3  (4) 


Whenever  f  =  tt/2  is  a  boundary  curve  for  the  integral  the  relation 
between  ft  and  Aamax  (j8,  /3n)  is  obtained  from  (2)  set  equal  zero,  giving 


cos  Aamax  (/3,  /3n)  =  -  ctn  j8jj  ctn /3 

Finally  we  obtain,,  when  £  =  n/2  is  a  boundary, 

*/3m  a 


(5) 


= 


/. 


cos  (8n  cos  j8  sin /3  cos-1  (—  ctn /3]\[  ctn  /3 ) 

dj8 


+  sin  /3n  sin2  /3  V  1  -  ctn2  /3 pj  ctn  /3 
While,  when  A  a  may  run  its  full  range  [0,  n  ]  (sin  A  a  =  0), 

Xfi  max  jfi  max 

cos  /3n  cos  /3  sin  /3  d  /3  =  cos  /3pj  sin2  /3 

Pm  in  P min 

To  integrate  (6)  we  first  simplify:  Let  A  =  sin  (3^  ,  B  =  cos  /3n 

Pn 

i 


(6) 


(7) 


Xr  m  a  x 

B  cos-1  ( - &-  Ctn  p\  d  (sin2  )8  ) 

Pm  in.  ' 


oPn 


~zr  f  V  A2  -  cos2  /3  d(cos  B) 


(8) 
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Integrating  by  parts  the  first  part  of  (8) 

I  fim  ax 


B 


-1  ctn  /3^  sin2  /3 


-  —  ^  B  sin2  /3  cl  |cos " 1  ctn  flj  J  -  ~~  J  V  A2  -  cos2  /3  d  ( cos  /3 )  (9) 


which  reduces  to 

=  “ 

S  77 


B  cos  1  | 

(  B 

—  -  ctn  p 

' 'j  sin2  6 

V  A 

I 

B2  d  (  cos  /3  ) 
V  A2  -  cos2  f3 


AfV~ 


A2  -  cos2  j3  d  ( cos  /3  )  (10) 

The  last  two  integrals  are  standard  forms  giving 


I 


B  cos  “ 1  — —  ctn  /3  ^  sin2  /3 


/3n 


B2  f  cos  (3 

-  sin"1  ( - 7 

77  \  A 


fim  a  x 


/3„ 


—  ^cos  /3  V  A2  -  cos2  /3  +  A2  sin"1  (  c 


By  definition  (prior  to  (8))  A2  +  B2  =  1 ,  hence 


U  =  ~ 

s  77 


B 


( - ctn  ft')  sin2  /3 


—  —  cos  |8  V  A2  -  cos2  /3 


/^m  ax 
/^m  in 


(ID 


and  this  is  defined  only  for  n/2  -  y  <  /3n  <  77/2  +  y  ,  /3  <  77/2  -  /Sn  • 
From  (11)  and  (7),  all  cases  are  derived.  The  integrations  are  abbre¬ 
viated  here.  (See  Figs.  Ill  - 1  and  III  -  2 .  ) 

Case  1.  I77  is  the  type:  It  s  0 


Case  2.  I2 


y  Aa(/3,/3N) 

-  2-  f  f 

**  fimin  ~  Z^N  77/2  «/o 


I^r  is  the  type. 
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Case  3. 

pfimin  =  ^/2  ~/3N  r>rr 

U  -  v  /  J 

77  "o  *'o 

_Aa(/3,/3N) 

—  +  i  L  f 

Pm  in  0 

(of  type  ln ) 

(of  type  l£  ) 

Case  4. 

i  y  _  n 

I4  =  -| r  J  J  -  (of  type  Irr  ). 
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•  y 


Case  1.  pN>  J+  Y 

None  of  Horizon  Circle  in  View  from  N. 

Il  a  0  for  7T  >  PN>  y  +  y 
Case  1,  2  I  ^  2  =  0  P|\|  =  ~2  +  ^ 


p=y 

„  7T 
3=2 


B=y 


7T 


Case  2.  N  in  Hemisphere  of  6  =  7T 
rr  tt 

2  5  3n  <  2  /  y 


1 


\2  =  7  cos  PN  sin^y 


/  .  -ifi 

sm  51 


sin^BN-cos^y 


sinBNcosy 


1  [sin-1 501L 
tt  \  sin  pN 


+  cos  y  ^si  n2B  N-cos2y] 


Case  2,  3 

I 


0 


TT 


2,3 


TT 


N  "  2 
y  -  sin  y  cos  y 


/ 


Case  3.  N  in  Hemisphere  of  6  ■  0 


tt  ^  „  tt 

2  5  PnI  2 


1  o  (  Jsin\-cos2r\ 

u  cos  PNsin^y  I  tt-  sin  1  ¥-r— * - ) 


sin  pN  cos  y  j 


1  /  .  -i  cos  y  r7T„  T 

H  \s,n  sTnpTT  +  cos  y  Vsm  Pn"cos  V  1  + 


TT 


Case  3,  4  PN  -  ~f  ~  Y 
l3j  4  ■  sin3  y 

Case  4.  All  Horizon  Circle  in  View  from  N. 


TT 


o<Pn<T"Y 


1 4  ■  cos  PN  sin^y 


Case5.  Pn"0 


sinT 


ig.  1 1 1-1  Cases  of  the  Earth  Radiance  integral 


Fig.  111-2  Horizon  and  Normal  Limits  on 
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APPENDIX  IV 

SAMPLE  PROFILE 


To  illustrate  the  nature  of  the  albedo -earth  radiance  profiles,  we 
take  as  model  a  cone-cylinder  vehicle,  in  a  circular  orbit,  tangent 
oriented.  The  cone  half-angle  is  taken  to  be  27  deg,  and  the  orbit  alti¬ 
tude  is  100  miles  above  the  earth  surface,  assumed  4000-mile  radius. 
The  orbit  is  inclined  33  deg  to  the  sun  direction  (Fig.  IV- 1). 

We  plot  the  earth  radiance  profile  integrals  in  Fig.  IV- 2  for  cone 
and  cylinder.  These  are  constant  functions,  independent  of  time,  for 
the  orbit  model  chosen. 

In  Fig.  IV- 3,  the  cylinder  albedo  profiles  are  shown  for  <£v  =  o  , 

90  deg.  By  symmetry,  these  functions  repeat  in  reverse  from  90  to 
180  deg.  Note  that  the  function  tails  vanish  at  about  ±  ( w/2  -  y)  of 

=  0  (j8n  =  77);  further  the  functions  become  everywhere  zero  for  about 
±y  of  0V  =  270  deg,  i.  e.  while  eclipse  occurs. 

In  Figs.  IV-4,  IV-5,  the  cone  albedo  profiles  are  shown  for^v  =  0, 
180  deg.  Because  the  receiving  surface  is  a  cone,  symmetry  fails,  and 
even  the  top  of  the  c  ine  receives  illumination  (the  tails  do  not  vanish). 

These  functions  were  obtained  by  numerical  integration  as  men¬ 
tioned  in  section  4.  3. 
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a.  Orbit  Model 


Dimensions  of  Apollo-type 
Vehicle  (Fig.  3) 


b.  Vehicle  Model 
Fig.  IV-1 
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Cone  Cylinder 


Fig.  IV-2  Earth  Radiance  Profiles 
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Fig.  IV -3  Albedo  —  Cylinder  Profiles 
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